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ABSTRACT 

We develop and test a model for the cosmological role of mergers in the formation and quenching of red, 
early-type galaxies. By combining theoretically well-constrained halo and subhalo mass functions as a function 
of redshift and environment with empirical halo occupation models, we predict the distribution of mergers as a 
function of redshift, environment, and physical galaxy properties. Making the simple ansatz that star formation 
is quenched after a gas-rich, spheroid-forming major merger, we demonstrate that this naturally predicts the 
turnover in the efficiency of star formation and baryon fractions in galaxies at ~ (without any parameters 
tuned to this value), as well as the observed mass functions and mass density of red galaxies as a function 
of redshift, the formation times of early-type galaxies as a function of mass, and the fraction of quenched 
galaxies as a function of galaxy and halo mass, environment, and redshift. Comparing to a variety of semi- 
analytic models in which quenching is primarily driven by halo mass considerations or secular/disk instabilities, 
we demonstrate that our model makes unique and robust qualitative predictions for a number of observables, 
including the bivariate red fraction as a function of galaxy and halo mass, the density of passive galaxies at high 
redshifts, the emergence/evolution of the color-morphology-density relations at high redshift, and the fraction 
of disky/boxy (or cusp/core) spheroids as a function of mass. In each case, the observations favor a model in 
which some mechanism quenches future star formation after a major merger builds a massive spheroid. Models 
where quenching is dominated by a halo mass threshold fail to match the behavior of the bivariate red fractions, 
predict too low a density of passive galaxies at high redshift, and overpredict by an order of magnitude the mass 
of the transition from disky to boxy ellipticals. Models driven by secular disk instabilities also qualitatively 
disagree with the bivariate red fractions, fail to predict the observed evolution in the color-density relations, and 
predict order-of-magnitude incorrect distributions of kinematic types in early-type galaxies. We make specific 
predictions for how future observations, for example quantifying the red fraction as a function of galaxy mass, 
halo mass, environment, or redshift, can break the degeneracies between a number of different assumptions 
adopted in present galaxy formation models. We discuss a variety of physical possibilities for this quenching, 
and propose a mixed scenario in which traditional quenching in hot, quasi-static massive halos is supplemented 
by the strong shocks and feedback energy input associated with a major merger (e.g. tidal shocks, starburst- 
driven winds, and quasar feedback), which temporarily suppress cooling and establish the conditions of a 
dynamically hot halo in the central regions of the host, even in low mass halos (below the traditional threshold 
for accretion shocks). 

Subject headings: quasars: general — galaxies: active — galaxies: evolution — cosmology: theory 



1. INTRODUCTION 

Recent, large galaxy surveys such as SDSS, 2dFGRS, 
COMBO- 17, and DEEP have demonstrated that the lo- 
cal distribution of galaxies is bimodal with respect to a 
number of physical properties, including color, morphol- 
ogy, star formation, con centration, and surface brightness, 
(e.g. IStrateva et al.l 1200 Ih . and that this bimodality extends 
at least to moderate redshifts, z ~ L5 (e.g.. Bell et al. 2004; 
IWillmer et all |2006|) with a significant population of mas- 
sive, red, passivel y evolving galaxies at even higher red- 
shifts (iLabbe et al.l l2b05: Krie k et al.ll2006h . The massive red 
galaxies in this bimodal distribution correspond to traditional 
spheroids, with high s urface brightness and concentration 
(iKauffmann et al. 2003). with little continuing star formation 
since their formation at early times (Trager et al. 2000). Un- 
derstanding the formation, and in particular the turning off or 
"quenching" of star formation on the red sequence, is there- 
fore of fundamental importance to understanding the origin of 
galaxies. 

Hierarchical theories of galaxy formation and evolution 
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indicate that large systems are built up over time through 
the merger of smaller progenitors, and galaxy interactions 
in the local Universe motivate the " merger hypothesis" 
jToomre & Toomr^l 19721: lToomreiri977l) , according to which 
colUsions between spiral galaxies produce the massive ellip- 
ticals observed at present times. 

Observations increasingly support the notion that galaxy 
mergers produce starbursts and structure ellipticals. The most 
intense starbursts, ultraluminous infrared galaxies (ULIRGs), 
are al ways associated with mergers (e.g. ISanders & Mirabell 
1996), with dense gas in their centers providing material to 
feed black hole (BH) growth and to boost the concentra- 
tion and central phase space density to match those of el- 
lipticals dHernquist et al.lll993l: [Robertson et al.ll2006l) . Like- 
wise, observations of in dividual merging s^steni s and gas-rich 
merger remnants (e.g.. Lake & Dressier 1986; Dovo n et al.l 
1994; Shier & Fischer 1998; James et al. 1999), as well as 
post-starburst (Eh-A/Kh-A) galaxies (Goto 2005), have shown 
that their kinematic and photometric properties, including ve- 
locity dispersions, concentrations, stellar masses, light pro- 
files, and phase space densities, are consistent with their 
eventual evolution into typical ^ elliptical galaxies. The 
correlations obeyed by these mergers and remnants (e.g.. 
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iGenzel et all 120011: iRothberg & JosephI l2006allbl) are similar 
to e.g. the observed fundamental plane and Kormendy rela- 
tions for relaxed ellipticals, and consistent with evolution onto 
these relations as their stellar populations age. This is fur- 
ther supported by the ubiquitous presence of fine structures 
such as shells, ripples, and tidal plumes i n ellipticals (e.g. 
ISchweizer & Seitzedl992l:ISchweizedl996h. which are signa- 
tures of mer gers (e.g. 'Ouinnl ll984tlHernquist & Ouinnll 19871; 
lHernquist & Spergel 199^rand the clustering and mass den- 
sity of ellipticals, consistent with passive evolution after for- 
mation in mergers (Hopkins et al. 2007d|) . 

Numerical simulations performed during the past twenty 
years verify that major mergers of gas-rich disk galaxies can 
plausibly account for these phenomena and elucidate the un- 
derlying physics. In iHopkins et al.l (l2007b). we provide an 
outline of the phases of evolution that might be associated 
with a major merger in the lifetime of a massive galaxy, but 
we briefly summarize them here. Tidal torques excited dur- 
ing a me rger lead to rapid inflows of gas into the centers o f 
galaxies (lHernquistl ll989; Barnes & Hernquist 1991,1991), 
triggering starbursts (Mihos & Hernquist 1994b, 1996) and 
feeding rapid black hole growth jPi Matteo e t al. 2005). Gas 
consumption by the starburst and dispersal of residual gas by 
supernova-driven win ds and feedback from black hole growth 
dSpringel et al.l l2005a) terminate star formation so that the 
remnant quickly evolves from a blue to a red galaxy. The stel- 
lar component of the progenitors provide s the bulk of the ma- 
terial for produc i ng the remnant spheroid jBarnesI 198811 19921 : 
lHernquislll992i Il993h through violent relaxation. A major 
merger is generally required in order for the tidal forces to ex- 
cite a sufficiently strong response to set up nuclear inflows 
of gas and build massive spheroids. Although simulations 
suggest that the precise definition of a major merger in this 
context is somewhat blurred by the degeneracy between the 
mass ratio of the progenitors and the orbit of the interaction 
jHernquist 1989; Hernquist & Mihos 1995; Bournaud et al. 
120051). svstematic st udies with both nume rical simulations 
dYounger et al.l |2007|) and observations (.Dasvra et al.l 120061 ; 
IWoods et al.ll2006l) " find that strong gas inflows and morpho- 
logical transformation are typically only observed below mass 
ratios 3 : 1, despite the greater frequency of higher mass- 
ratio mergers. In what follows, unless explicitly noted, we 
generally mean the term "mergers" to refer specifically to ma- 
jor mergers. 

It also must be emphasized that essentially all numer- 
ical studies of spheroid kinematics find that only merg- 
ers can reproduce the observed kinematic proper ties of ob- 
served elliptical galaxies and "classical" bulges (iHernquistI 
1989, 1992, 1993; Barnes 1988, 1992; Schw eized I1992F 
Naab & Burkert_j2003t Bournaud et al. 2005; Naab et alj 
2p06a b; Naab & TruiiUd 2006; Jesseit et al. 2006; Cox et al] 
2006b). Disk instabilities and secular evolution (e.g. bar 
instabilities, harassment, and other isolated modes) can in- 
deed produce bu l ges, but these are invariab l y "pseudobulges" 
dSchwarzl [T9p; lAtha nassoul a et all 1198 3'; 'Pfennigei"1984t 
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lAthanassou la 2005), with 
clearly distinct shapes (e.g. flattened or "peanuf'-shaped 
isophotes), rotation properties (large v/cr), internal correla- 
tions (obeying different Kormendy and Faber-Jackson rela- 
tions), light profiles (nearly exponential Sersic profiles), and 
col ors and/or substructure from classical bulges (for a review, 
seelKorm endv & Kennicutt 2004) . Observations indicate that 
pseudobulges constitute only a small fraction of the total 



mass density iri spheroids (< 10% ; see lAUen et aTl l2006al: 
iBaUet al.1 120061; iDriver et al.H2007h . becoming a large frac- 
tion of the bulge population only in small bulges in late-type 
ho sts (e.g. Sb/c, corresponding to typical Mgai 5 . IO^^Mq; 
see lCarollo et al.ll998HKormendy & Ken nicutt2004. and ref- 
erences therein). This is not to say that secular processes 
cannot, in principle, build some massive bulges (see e.g. 
iDebattista et al.ll2004 1200 6"). However, although such pro- 
cesses may be important for the buildup of low mass black 
hole and spheroid populations, it is empirically clear that sec- 
ular evolution cannot be the agent responsible for the forma- 
tion of most elliptical galaxies. 

Motivated by these considerations. iHopkins et aP (I2006all3) 
developed a model where starbursts, quasars, supermassive 
black hole growth, and the formation of red, elliptical galax- 
ies are connected through an evolutionary sequence, caused 
by mergers between gas-rich galaxies. It is important to keep 
in mind that this does not rule out other processes occurring at 
lower levels and under other circumstances. For example, we 
are not claiming that all bulges result from mergers - secu- 
lar pseudobulge growth does appear to be important for small 
bulges in disk-dominated systems, and additional processes 
may act to redden satellite galaxies in massive halos, a poten- 
tially important contributor to the population of red g alaxies at 
low masses Mgai < lO'^M© (e.g. lBlanton et alj|2005b . More- 
over, spheroid evolution by gas-free ("dry") mergers will go 
on, but does not explain how stellar mass is initially moved 
onto the red sequence or transformed from disk to spheroid. 

All of this, however, only goes to the question of the 
formation of elliptical galaxies, not to the question of how 
such galaxies become (and stay) "red and dead." It is well 
establi shed from both numer ical simulat ions (Springe l et alj 
l2005ah and observations (e.g. lRothberg & Joseph 2006ah that 
merger remnants redden rapidly onto the red sequence as typ- 
ical early-type galaxies. However, it is still debated whether 
or not such systems will stay on the red sequence for long 
periods of time, since this requires some suppression of sub- 
sequent accretion and cooling as their host dark matter halos 
grow. In massive elliptical galaxies, it is not obvious how the 
formation of cooling flows has been suppressed since z ^ 2, 
despite observations finding that the cooling times of large 
quantities of gas are shorter than a Hubble time. In other 
words, there is an important outstanding question, which we 
seek to address: do major mergers or their remnants effec- 
tively quench future star formation (i.e. maintain low star for- 
mation rates for significant cosmic times), or is it some other, 
independent process which is responsible for quenching? 

At low redshift, there appears to be a clear association 
between quenched (red, passive) galaxies and the presence 
of a massive spheroid, at least for the relativel y massive 
Mg3L> 10^ °Mfn systems o f interest in this paper. iBell et afl 



( 120031) andlMcIntoslLet al. (2005) find that > 80% of the z = 
red population are classical, bulge-dominate d systems, with 
most o f the remainder being early-type disks. iDrory & Fished 
(l2007h further investigate these disk-dominated systems, and 
find that early-type disks on the red sequence have uniformly 
classical bulges (presumably formed via mergers), whereas 
disks of comparable mass, luminosity, and bulge size host- 
ing pseudobulges (formed via secular instabilities) remain in 
the blue cloud. At higher redshifts, morphological signatures 
are less clear, and an increasingly large fraction of red galax- 
ies (naively identified by simple color cuts) are contaminant 
dusty or edge-on disks (clearly not true quenched/passive sys- 
tems). However, those systems which can be clearly iden- 
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tified as truly passive appear to be overwhelmingly compact 
spheroids (Mcintosh et al. 2005; Bundv et al. 20( p), even at 
z~2-3 (Labbe et al., .2005.: .Kriek et al., .20061: Izirm et alj 
|2007|) . This suggests a strong connection between a major, 
spheroid-forming merger and galaxy quenching. 

The standard framework for understanding quenching fol- 
lows the cooling of gas in the galaxy host halo. From simple 
scaling arguments one can show that at low halo masses the 
cooling time will (in the absence of heating mechanisms) be 
shorter than the free-fall time of the gas, and accretion is only 
limited by the free-fall of newly accreted halo gas onto the 
central galaxy - the so-called "rapid cooling" or "cold accre- 
tion" regime. Once the halo becomes sufficiently massive, the 
cooling time becomes longer than the free fall time, and so gas 
does not simply fall onto the central galaxy, but rather forms a 
quasi-static, pressure supported hydrostatic equilibrium - the 
"hot halo" regime. New gas accreted will shock against this 
pressure-supported structure, heating itself and the gas inte- 
rior to it, and accretion will proceed only gra dually, from the 
cooling of the gas at the center of the halo ( Rees & Ostrikerl 
[1977; Norman & Silk 1979; Blumenthal et al. 198^^ 

Numerical simulations suggest t hat this transition occurs 
at a mass Mh.io ~ lO" - IO'^Mq (iBirnboim & Deke]|l2003t 
iKeres et akl l2005h . In many prescriptions (such as the 
"halo quenching" models to which we refer in § I3.2l l. it 
is simply assumed that the development of a hot halo at 
this mass threshold is the dominant criterion for quench- 
ing. How ever, both numerical simulatio ns and analytic cal- 
culations dkauffmann & Haehn elt 200rt iBenson et allllOOl 
iKeres et ani2005[ and references therein) argue that this tran- 
sition alone cannot solve the "cooling flow" problem - namely 
that the high densities at the core of the pressure-supported 
hot halo will allow rapid cooling onto the central galaxy, pro- 
ducing large galaxies which are much too massive, gas-rich, 
disk-dominated, actively star-forming, young and blue rela- 
tive to the observations. Some kind of heating term is needed 
to prevent this from occurring and maintain quenching. 

It has become popular to invoke activity from a low- 
Eddington ratio AGN in the central galaxy a s the source 
of this heating term - the "radio-mode" AGN (' Croton et al.l 
12006). This, however, requires the presence of a massive 
black hole (and therefore a correspondingly massive spheroid) 
accreting in a relatively low steady state (i.e. with most of the 
cold gas in the galaxy consumed). This requirement, along 
with the arguments above, suggests that merger history might 
be just as important as (if not more important than) halo mass 
in determining the quenching of a given galaxy. Ultimately, 
we emphasize that the detailed numerical simulations and an- 
alytic calculations of the hot halo regime do not argue that 
entering this regime does, or can, directly quench future cool- 
ing. Rather, these calculations argue only that the "hot halo" 
regime provides an ideal environment in which quenching 
mechanisms might operate. 

Unfortunately, obtaining a purely theoretical framework 
for any quenching scenario is difficult because cosmolog- 
ical simulations including gas dynamics currently lack the 
resolution to describe the small-scale physics associated 
with disk formation, galaxy mergers, star formation, and 
black hole growth. A popular alternative has been the 
employment of semi-analytic methods, adopting various 
prescriptions for quenching and feedback processes and 
comparing the predictions with observe d galaxy popula- 
tions ( e.g. iKauffmann & Haehneli 120001: ISomerville et al.l 
I2OOII IBenson et al.1 l2003t iKhochfar & BurkertI l2003t 
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models have robustly shown the need for some quenching 
processes, and their great success has been demonstrating 
that simple prescriptions for basic feedback elements yield 
good agreement with local g alaxy mass/lumino s ity functions 
and color distributions (e.g . ICroton et"a!] 120061 : iBower et al] 
l2006tlCattaneo et al.ll2006h . 

However, the similar success of a large variety of such 
prescriptions at matching these basic local constraints has 
demonstrated that such predictions are fundamentally non- 
unique. For example, ICattaneo et al.l (l2006h have shown 
that one obtains similar galaxy mass functions and color- 
magnitude relations whether one adopts a pure halo mass 
threshold for quenching, a halo mass threshold which de- 
pends on some feedback balance with a low-luminosity AGN, 
or a (halo mass-independent) galaxy bulge-to-disk criterion. 
Clearly, these simple constraints are insufficient to discrimi- 
nate between the mechanisms associated with galaxy quench- 
ing. Furthermore, the diversity of semi-analytic prescriptions 
has demonstrated that there are considerable degeneracies be- 
tween, for example, the prescriptions for star formation in 
disks and those for quenching, despite the fact that the two 
should be constrained by independent galaxy populations. It 
is therefore necessary to determine what, if any, are the robust 
differences between various quenching prescriptions, and to 
study higher-order observational constraints (such as e.g. the 
redshift evolution of populations, or bivariate distributions of 
galaxy properties as a function of both galaxy mass and halo 
mass or galaxy kinematics) that hold the potential to break 
these degeneracies. 

In th e first of a pair of companion papers dHopkins et alj 
I2007bl henceforth Paper I), we describe a strategy that en- 
ables us, for the first time, to provide a purely theoretical 
framework for our models of merger-induced activity. By 
combining previous estimates of the evolution of the halo 
mass function with halo occupation models and our estimates 
for merger timescales, we infer the statistics of mergers that 
form spheroids. Because our merger simulations relate star- 
bursts, quasars, and red galaxies as different phases of the 
same events, we can graft these simulations onto our theoreti- 
cal, cosmological calculation and determine the cosmological 
birthrate of these various populations and their evolution with 
redshift. In particular we demonstrate in what follows that 
there are a number of unique, robust predictions of a model 
in which mergers drive the quenching of galaxies (in addition 
to forming spheroids in the first place), distinct from the pre- 
dictions of models in which this quenching is set just by halo 
properties or secular (disk) instabilities. We find that obser- 
vations of red galaxies support our predictions, and disfavor 
other theoretical models. 

In Paper I, we describe our model and use it to investi- 
gate the properties of mergers and merger-driven quasar ac- 
tivity. In this paper (Paper II), we extend this to study the 
properties of merger remnants and the formation of the early- 
type galaxy population. We begin by briefly reviewing the 
key elements of the model from Paper I in § |2] In § [3] we 
use the method developed in Paper I to examine the conse- 
quences of a general model in which major merger remnants 
remain "quenched" once the merger terminates star forma- 
tion. Specifically, § 13.11 shows the predictions of this model 
for the buildup of early-type or red galaxy mass functions and 
mass density with redshift, and the formation times of early- 
type galaxies. In § l3.2l we demonstrate how the resulting frac- 
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tion of red or "quenched" galaxies depends on properties such 
as halo and galaxy mass, and contrast these with the predic- 
tions of alternative models in which the quenching is asso- 
ciated with a halo mass criterion or secular processes (disk 
instabilities). In § 13.31 we extend these comparisons to the 
redshift evolution of these trends. In § l3.4l we briefly examine 
the role of subsequent gas-poor major mergers in this model, 
and compare with observations of early-type galaxy structure. 
In §|4]we outline the broad physical mechanisms which give 
rise to such a model. We examine in § l4.1l how mergers are as- 
sociated with the "transition" of galaxies from the blue cloud 
to the red sequence, and in § l4.2l we examine the role of differ- 
ent feedback mechanisms in "maintaining" low star formation 
rates in remnant elliptical galaxies. We discuss and summa- 
rize our conclusions in § |5] 

Throughout, we adopt a WMAP3 (JIm, ^ h, a^, = 
(0.268, 0.732, 0.704, 0.776, 0.947) cosmology (ISoergel et all 
12006). and normalize all observations and models shown to 
this cosmology. Although the exact choice of cosmology may 
systematically shift the inferred bias and halo masses (pri- 
marily scaling with cts), our comparisons (i.e. relative biases) 
are for the most part unchanged, and repeating our calcu- 
lations for a "concordance" (0.3,0.7,0.7,0.9, 1.0) cosmol- 
ogy or the W MAPl (0.27,0.73,0.71,0.84,0.96) results of 
ISpergel et al.l (l2003i) has little effect on our conclusions. We 
also adopt a diet Salpeter IMF following Bell et al. (2003^), 
and convert all stellar masses and mass-to-light ratios to this 
choice. Again, the exact choice of IMF systematically shifts 
the normalization of stellar masses herein, but does not sub- 
stantially change our comparisons. UBV magnitudes are in 
the Vega system, and SDSS ugriz magnitudes are AB. 

2. MERGERS: THE BASIC MODEL 

The model which we use to calculate the rate and nature of 
mergers as a function of e.g. mass, redshift, and environment 
is described in detail in Paper I, but we briefly outline the key 
elements here. 

1. Halo Mass Functi on: We begin by a dopting the halo 
mass function following iSheth et al.l (1200 lb . There is little 
ambiguity in this cal culation at all redshifts and masses of 
interest (z < 6; e.g. Reed et al. 2007), and we do not consider 
it a significant source of uncertainty. 

2. Subhalo Mass Function: The subhalo mass function 
of each halo is then calculated. Although numerical simu- 
lations and semi-analytic calculations generally give similar 
results (especially for the major-merger mass ratios of inter- 
est in thi s paper, as opposed to ver y small subhalo popula- 
tions; see lvan den Bosch et al.ll2005h . there is still some (typi- 
cal factor < 2) disagreement between different estimates. We 
therefore repeat most of our calculations adopting b oth our 
"default" subhalo m ass function calculation (Zentn eF et alJ 
l200 5t iKravtso v et al.i 2004) and an alternative subhalo mass 
function calculation dvan den Bosch et af] 2005h (normaUze d 
to match cosmological simulations as in LShawet al.ll2006h . 
which bracket the range of a number of different estimates 
(e.g., Springel et al. 2001; Tormen et al. 2004; De Luci a et al] 
|2004; Gao et al. 2004; Nurmi et al. 2006) and demonstrate die 
uncertainty owing to this choice. The difference is ultimately 
negligible at Mgai > 10'"Mq (where, unless otherwise spec- 
ified, Mgai refers to the baryonic mass of the galaxy) at all 
redshifts, and rises to only a factor ^ 2 at Mgai < IO^'Mq 
(probably owing to differences in the numerical resolution of 
various estimates at low halo masses). 



3. Halo Occupation Model: We then populate the cen- 
tral galaxies and "major" subhalos with an empirical halo 
occupation model. Although such models are constrained, 
by definition, to reproduce the mean properties of the halos 
occupied by galaxies of a given mass/luminosity, there are 
known degeneracies between parameterizations that give rise 
to (typical factor ^ 2) differences between models. We there- 
fore again repeat all our calculations for our "default" model 
(IConroy et alJl2006l) (see also Vale & Ostriker 2006) and an 
alternate halo occupation model (Yang et al. 2003) (see also 
lYan et alj|2003l; IZheng et al. 2005 ), which bracket the range 
of a number of calculations (e.g., Ivan den Bosch et al.ll2006l: 
ICoorayll2005[ 120061: IZheng et alJl2005h . Again, we find this 
yields negligible differences atMgai > IO^^Mq (as the cluster- 
ing and abundances of massive galaxies are reasonably well- 
constrained, and most of these galaxies are central halo galax- 
ies), and even at low masses the typical discrepancy in our 
predictions owing to the choice of halo occupation model rises 
to only ^ 0.2 dex. 

We note that we have also considered a variety of pre- 
scriptions for the redshift evolution of the halo occupation 
model: including that directly prescribed by the quoted 
models, a complet e re- derivation of the HOD models of 
IConrovetal] ( l2006h and I Vale & Ostriked (l2006| ) at different 
redshi ft s from th e observed mass functions of Bundv et alJ 
(j20()5h : iFbntaiia etal. (2006); Borchetal. (2006); Blanto^ 
(2006) (see Paper I), or assuming no evolution (in terms of 
galaxy mass distributions at fixed halo mass; for either all 
galaxies or star-forming galaxies). We find that the resulting 
differences are small (at least at z < 3), comparable to those 
inherent in the choice of halo occupation model. This is not 
surprising, as a number of recent studies suggest that there 
is little evolution in halo occupation parameters (in t erms of 
mass, or relative to L^,) wit h redshift (lYan et alJl2003tlCoorayl 
120051; IConrov et al.l l2006h , or equivalently that the masses 
of galaxies hosted in a halo of a given mass are primarily 
a function of that halo ma ss, not of redshift (IHeymans et au 
l2006al: IConrov et ani2007l) . This appears to be especially true 
for star-forming and ^ L» galaxies (of greatest importance 
for our conclusions; 'Conroy e t alJl200'7h . unsurprising given 
that quenching is not strongly operating in those systems to 
change their mass-to-light ratios. 

4. Merger Timescale: Having populated a given halo and 
its subhalos with galaxies, we then calculate the timescale for 
mergers between major galaxy pairs. This is ultimately the 
largest source of uncertainty in our calculations, at all red- 
shifts and masses. Again, we emphasize that some of our 
calculations are completely independent of these timescales. 
However, where adopted, we illustrate this uncertainty by 
presenting all of our predictions for three estimates of the 
merger timescale: first, a simple dynamical friction for- 
mula (this is what is generally adopted in semi-analytic 
models, for example). Second, a group capture or colli - 
sional (i.e. effective gravitational) cross section (e.g. IWhitd 
1976"; 'Ki-ivitskv & Konto rovichI \l99% iMakino & Hull \l99% 
Mamon 2006) approximation, generally more appropriate on 
small scales, in satellite-satellite mergers, or in the merger 
of two small field halos. Third, an angular momentum (or- 
bital cross section) capture estimate (i.e. considering cap- 
ture into the effectiv e angular- momentum space of mergers; 
iBinnev & Tremaindfr 987). 

At large masses and redshifts z < 2.5, this is a surpris- 
ingly weak source of uncertainty, but the estimated merger 
rates/timescales can be different at low masses M„ai ^ 
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10 Mq and the highest redshifts z ^ 3-6. At low masses, 
this owes to a variety of effects, including the substantial dif- 
ference between infall or merger timescales and the timescale 
for morphological disturbances to be excited (different in e.g. 
an impact approximation as opposed to the circular orbit de- 
cay assumed by dynamical friction). Note that where relevant, 
we have used numerical simulations to estimate the typical 
duration of the final merger stages or e.g. the morphological 
relaxation time (in which mergers will be id entified by typica l 
morphological classification schemes, see Lotz et alj|20()7h . 
The difference in redshift evolution is easily understood: at 
fixed mass ratio, the dynamical friction timescale scales as 
fdf oc fn oc p"'/^, but a "capture" timescale will scale with fixed 
cross section as f oc l/(n (cr v)) cx p"' , so that (while the details 
of the cross-sections make the difference not quite as extreme 
as this simple scaling) the very high densities at high redshift 
make collisional merging grow rapidly in efficiency. The true 
solution is probably some effective combination of these two 
estimates, and the "more appropriate" approximation depends 
largely on the initial orbital parameters of the subhalos. At 
present, we therefore must recognize this as an inherent un- 
certainty, but one that serves to bracket the likely range of 
possibilities at high redshifts. 

In Paper I (§ 2.2), we show that together, these criteria nat- 
urally define a preferred major-merger scale (host halo mass 
Mhaio) for galaxies of mass Mgai - the "small group scale," 
only slightly larger than the average halo hosting a galaxy of 
mass Mgai. This is the scale at which the probability to accrete 
a second galaxy of comparable mass ~ Mgai (fuel for a major 
merger) first becomes significant. At smaller (relative) halo 
masses, the probability that the halo hosts a galaxy as large 
as Mgai declines rapidly. At larger masses, the probability that 
the halo will merge with or accrete another halo hosting a 
comparable ^ Mgai galaxy increases, but the efficiency of the 
merger of these galaxies declines rapidly. We stress that this 
small group scale is distinct from the more typical large group 
scale identified observationally (the average small group halo 
will still host only 1 galaxy of mass ^ Mgai, and groups will 
only consist of 2-3 members of similar mass). This is not 
to say, however, that mergers occur (in a global sense) at a 
specific scale, since the small group scale is different for dif- 
ferent galaxy masses - a consequence of this model is the ob- 
servational fact that mergers occur in halos of all masses and 
in aU environments (including field and even void environ - 
ments; LSol Alonso et ani2006li lGotoll20()l iHogg etal\f2(M) . 
although the characteristic masses and star formation histories 
of galaxies merging will change in different environments. 

In Paper I we compare this model with a number of obser- 
vations, and show that it reproduces the mass functions and 
star formation histories of galaxies, merger mass functions 
(and infrared luminosity functions) and merger fractions as 
a function of galaxy and/or halo mass and redshift, the clus- 
tering of mergers as a function of mass and redshift, and the 
dependence of merger rates and fractions on small-scale envi- 
ronmental properties. This provides some reassurance that we 
are accurately predicting the rate and nature of major mergers 
as a function of these properties, and can use this model to 
make robust predictions for the nature of merger remnants. 

3. ELLIPTICALS 

We now turn to the possibility of an association between 
mergers and the termination or quenching of star formation in 
remnant galaxies. In §|4]we consider potential physical mech- 
anisms for this quenching, but we caution that at present these 
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Fig. 1. — Predicted local quenched/red/early-type galaxy mass function 
(lines) obtained by integrating forward the major merger mass function to 
z = (i.e. assuming that each merger leaves a quenched early-type remnant). 
Different styles show different variants of our calculation which bracket the 
range of our uncertainties, varying e.g. the subhalo mass functions, halo oc- 
cupation model, and approximation used to calculate merger timescales (as 
described in § |2] see Paper I for a more detailed comparison). We com- 
gare with observed early-type or red galaxy mass functions from lBell et al.l 
(2003, SDSS elliptical and red galaxy mass functions; black and red X 's, re- 
spectively). Wake et al. (2006, SDSS LRGs; blue squares), and Jones et aD 
(2006, 6dF LRGs; pu rple diamonds). The mass functions from Wake et ay 
(2006) and I Jones et al. (2006) are converted from luminosity functions us- 
ing the luminosity-dependent mass-to-light ratios from lBell et al.it200 3). We 
show both the directly measured Wake et al. (2003) result (open) and that 
corrected for passive evolution from z = 0.i (filled). 

mechanisms are neither well-understood nor observationally 
well-constrained. As a consequence, we first wish to examine 
the consequences of the simple hypothesis that some mecha- 
nism quenches star formation after a major merger, whether 
it involves gas exhaustion, starburst or quasar feedback, hot 
halo formation, or other mechanisms. We therefore make the 
simple ansatz: Systems are quenched after a major merger of 
sta r-fo rming/gas-rich ga laxies . 

3.1. Integrated Populations 

In Paper I we calculated the major merger rate of galaxies 
as a function of galaxy mass and redshift. If each such merger 
leaves a quenched early-type remnant, then we can integrate 
the merger rate forward in time to obtain the early-type or red 
galaxy mass function at each redshift, 

r , df 

(/>early(Mgal) = / n(Mgal I ^ ^Z. (1) 

Figures [T]&|2] show this at several redshifts for our model of 
major mergers. Note that Equation [T] adds the contribution 
from all mergers - i.e. implicitly includes in the mass func- 
tion the contribution from "dry" or spheroid-spheroid merg- 
ers. Technically, we should also include the sink term from 
dry mergers, -2 /ndry(0.5Mgai |z)df, representing the loss of 
two early-types of mass ^ Mgai/2 for each major dry merger 
of final mass Mgai. This requires a number of additional as- 
sumptions for the red/blue galaxy fraction as a function of 
Mgai or Mhaio and the initial mass ratios of mergers, so we 
have not included it here, but note that for reasonable empir- 
ical estimates (such as those in § 13.21 ) of these numbers, the 
sink term has little effect. That is not to say that at low red- 
shift, the dry merger contribution cannot indeed be important 
to the shape of the mass function where it is falling steeply 
at high mass (;» M,). Because of this steep fall-off, moving 
a small fraction of lower mass systems to higher masses can 
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Fig. 2. — As Figure [T] but at each of several redshifts. Points at z = 
are as in Figure \l\ At h i gher redsh ifts, early-type or red galaxy MFs are 
shown from Bundy et al. 1 2005, 2006, red circles), Borch et al. 1 2006, purple 
squares), Franceschini et al. (2006, cyan stars), Pannella et al. (2006, orange 
triangles), Fontana et aU 12004 . blue inverted triangles). Wake et alj t2006t 
blue filled squares), and points at z = 1.7, 2.5, and 3.1 are estimated from the 
number density of passively evolving (non-star forming) red galaxies with 
stellar masses > 1O"M0 injDaddi et al. (2005, cyan square), ILabbe et al.l 
f200f, green square), and IGrazian et aU ( 120 07, magenta star), respectively. 
Masses (or mass ranges) have been corrected to our adopted IMF. The inte- 
grated merger mass function is consistent with the observed red galaxy mass 
function at all redshifts z ^ 0-3.5. 

significantly increase the number density of the most massive 
systems. However, the loss of less massive systems is a small 
correction. The dominant term at masses (< a few M,) impor- 
tant for the total mass density of red systems is the movement 
of systems to the red sequence by gas-rich mergers. 

We have also neglected growth via minor mergers: how- 
ever, we demonstrate in Paper I that this is also a small correc- 
tion; i.e. mass growth is dominated by major merg ers and star 
formation, as seen in cos mological simu lations ( Mailer et alj 
I2OO6,) and observations ( Zheng et al.l2007.) (although it is pos- 
sible that minor mergers become important for the most ex- 
treme, massive BCGs). 

As discussed in § |2] there are a number of uncertainties 
at the lowest masses Mgai < 10"' Mq, which are evident in 
the differences between our predictions in Figures [T] & |2] - 
these include issues of completeness and resolution in the 
subhalo mass functions and halo occupation models, and sen- 
sitivity (for very low-mass mergers) to the method used to 
calculate merger cross sections (for example, the difference 
between a dynamical friction and an impact approximation 
becomes large). The predictions in this regime are probably 
subject to a number of other caveats, as well. At the low- 
est masses Mgai < a few x 10^ M©, satellite-satellite mergers 
(the dynamics of which are sensitive to orbital parameters) 
become an important contributor to the total merger remnant 
population. Also, the fraction of observed pseudobulges starts 
to become large, implying that secular instabilities may be- 
gin contributing significantly to the early-type population be- 
low these masses. Finally, many of the observed red galax- 
ies at masses below this threshold (almost an order of mag- 
nitude below M*) are satellites of more massive systems, so 
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Fig. 3. — Integrated stellar mass density of red or early-type (passive) 
galaxies as a function of redshift, compared to the predicted stellar mass den- 
sity which has undergone a major galaxy merger (i.e. integrating the merger 
rate function of star-forming galaxies). Points are as in Figures [T]&|2] with 
lines showing the predictions from our different methods of calculating the 
merger rate function from §|2](see Paper I). We add the mass densi ty esti- 
mates from I Abraham et al . (2007, violet diamonds) and van Dokk um et all 
120061 orange square). Dotted green (uppermost) line shows the total stel- 
lar mass density of the universe expected from the integrated star formation 
history in Hopkins & Beacom (2006) (normalized to the z = value from 
IBell et al.l t2003i) ). The mass density of systems which have undergone ma- 
jor, gas-rich mergers agrees well at all redshifts with the mass density in red 
or early-type galaxies, with sufficient mergers occurring at high red shifts to 
account for the observed densities from Labbe et al. ( 2005) and Gra zian et al.l 
<2007l) atz~2-4. 

processes like ram pressure stripping, tidal stripping, harass- 
ment, and a cutoff of new accretion are likely to be impor- 
tant (and may even dominate their becoming red in the first 
place). At higher masses Mgai > IO'^Mq, however, the agree- 
ment between our predictions is good, regardless of which 
subhalo mass functions, halo occupation models, or merger 
timescale approximations we adopt. Moreover, almost all of 
these galaxies are obse rved to be central halo galaxies (e.g. 
IWeinmann et al.l2006al) and the pseudobulge fraction is small, 
so we can have some confidence that satellite and secular pro- 
cesses are not a large effect (see also § 13. 2b . We refer to Paper 
I for a more detailed comparison, but we have tested the model 
extensively for these masses (Mgai ^ IO'^'Mq), and find it is 
both robust and consistent with observed statistics of mergers 
as a function of galaxy and halo mass, redshift, and galaxy 
color/morphological type. 

Figure |3] plots the integrated version of this, namely the 
mass density of red/early-type systems as a function of red- 
shift. Here, we integrate only the gas-rich merger rate func- 
tion (the same merger rate function we used to predict the 
quasar luminosity function in Paper I; i.e. using our empiri- 
cal halo occupation model to identify specifically mergers of 
gas-rich or star-forming galaxies), as dry mergers cannot, by 
definition, increase the mass density of red galaxies. 

The integrated mass which has undergone major, gas-rich 
mergers agrees well with the mass density of red galax- 
ies at all redshifts. Even at high redshifts z ~ 2-4, this 
merger-driven model has no difficulty accounting for the rel- 
atively large mass densities of red galaxies observed by e.g. 
Labbe et alJ (12005), van Dokk um et al.l (l2006l) . iGrazian et al] 
(i2007l) . and iKriek et al.. (.200C as the highest-overdensity 
peaks in the early universe undergo rapid major mergers (sug- 
gested in the observations as w ell, given the color-density re- 
lation of these objects; Ouadri et al.ll2007h . It is important 
to note that a significant number of these high-redshift sys- 
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Fig. 4. — Predicted ages of early-type galaxies (at z = 0) as a function of 
stellar mass from the integrated mass functions in Figure[2] compared to ob- 
servations from Nelan et al. (2005, red circles), Thomas et al. ( 2005, orange 
squares; we take their mean values as opposed to those in a specific environ- 
ment), and Gallazzi et al. ( 2006, blue stars). Errors show the mass ranges and 
dispersion in ages within each mass range (not the error in the mean ages). 
Blue lines show the predicted mean lookback time to the final gas-rich merger 
for different estimates of the merger rate as in Figure[T] Black Hnes show the 
same, but with age calculated from the mean redshift (as opposed to lookback 
time), showing the systematic offset owing to this choice of definition. 

tems have been spectroscopically con firmed as passive , "red 
and dead" systems (Ki-iek et al. 2006; W uvts et aO2007h with 
elliptical morphologies (Zirm et aL .2007 and relatively old 
ages a few x 10^ yr, or ~ 1 /Sfe at these redshifts). This 
is in strong contrast to some pure hot-halo quenching mod- 
els, in which cold accretion within a hot halo persists at high 
redshifts z > 2. 

We consider a detailed comparison wi th these models in 
§ 13.31 but note, for example, that the iDekel & BirnboirnI 
(|2006|) estimate of the hot halo quenching mass predicts that 
at z 3.5, cold flows continue within all Mhaio ^ 10'"* M0 ha- 
los, which allows only a completely negligible maximum red 
galaxy mass density (even if we adopt a 100% baryon con- 
version efficiency and assume all quenched halos are "red") 
- even lowering this threshold by an order of magnitude pre- 
dicts a quenched galaxy mass density an order of magnitude 
below that observed. 

Having integrated forward the implied rate of formation of 
early-type galaxies, we can also predict the ages of early types 
as a function of their mass. Figure |4] shows this compari- 
son. We note that age here is defined as the time since the 
last gas-rich major merger (systems may, of course, undergo 
subsequent gas-free mergers, but this will not contribute new 
star formation). Our model contains no information about the 
prior star formation histories of merging disks. However, the 
observations to which we compare typically measure single 
stellar population (single burst) or light-weighted ages, which 
tend to reflect the last significant epoch of star formation. 

We emphasize that this does not imply that most of the stars 
in spheroids form in a short-lived, merger-induced burst. Di- 
rect calculation of the inferred stellar po pulation ages from 
line index and SED fitting (following Trager et al.' 2000) for 
realis tic star formation histories from the semi-analytic mod- 
els of lSomerville et all (1200 Ih and the h ydrodynamical merger 
simulations of [Robertson et al.l (1200 6') suggests that the ages 
inferred for present early-type galaxies indeed reflect the 
epoch of the termination of star formation, even when > 95% 
of stars are formed over a much longer timescale at signifi- 
cantly earlier times (in these cases, in quiescent star formation 



in disks). 

Since we are interested in testing the possibility that major, 
gas-rich mergers are associated with the termination of star 
formation, these ages are the most appropriate with which to 
compare. But we again emphasize the caveat that without 
the details of the star formation histories in progenitor disks, 
our ages are subject to some systematic uncertainties. In any 
case, the agreement is good, suggesting that mergers have the 
correct timing, in a cosmic sense, to explain the shutdown of 
star formation in early-type systems. 

3.2. Color-Density Relations: The Dependence of Red 
Fractions on Halo Mass and Environment 

We now study the distribution of red galaxies in greater 
detail, to highlight the unique features of a merger-driven 
quenching model. For clarity, we focus only on central galax- 
ies, and ignore the (potentially) completely physically distinct 
mechanisms (ram pressure stripping, tidal heating, etc.) re- 
sponsible for quenching satellite galaxies in massive halos. 
Therefore, in what follows, our comparison of quenching and 
red/blue galaxies explicitly ignores satellite galaxies. 

Figure |5] illustrates three qualitatively distinct classes of 
models for quenching. We distinguish our "merger" model 
(systems quench after a major, gas-rich merger), a pure "halo 
quenching" model (systems quench upon crossing a critical 
halo mass), and a "secular" model (internal processes - set 
by the galaxy mass and/or size - solely determine galaxy 
color/star formation history). The models all predict that the 
most massive halos and most massive galaxies are predomi- 
nantly quenched. However, in detail, the models differ in the 
behavior of quenching with respect to galaxy and halo mass. 

In the simplest halo quenching models, the ability of a 
galaxy to redden is completely determined by its host halo 
mass. In the simplest secular models, this is completely de- 
termined by the galaxy mass. In contrast, a merger-driven 
model depends on both - mergers will proceed more rapidly 
and efficiently at high Mgai in a given Mhaio, and larger Mhaio 
systems represent larger overdensity peaks which are more 
evolved and more likely to have undergone a period of merg- 
ing (recall, we refer to accreting a pair of mass ^ Mgai to fuel 
a major galaxy merger). 

Of course, the cartoon illustration in Figure |5] ignores some 
details. In many halo quenching models, quenching also re- 
quires a massive BH or some other feedback mechanism, 
which implicitly requires a relatively massive spheroid and 
therefore depends to some extent on the stellar mass and 
merger history of the system. In many secular models, galaxy 
structure and disk instability are influenced by halo properties 
(e.g. concentration) that vary with halo mass and accretion 
history. We therefore consider a more detailed comparison 
with sta te-of-the-art semi-an alytic models. We extract the re- 
sults of lCroton etaP (l2006l) ^ and iBower et"an (l2006l) l both 
recent fully cosmological semi-analytic models based on the 
Mille nium dark-m atter simulation (Springel et al. 2005c). 

The lCroton et alJ (120061) models correspond roughly to the 
halo quenching models described above - a massive BH 
is required to maintain the hot halo, but development of 
the hot halo reservoir (upon crossing the appropriate halo 
mass threshold) is still effec tively the dominant criterion 
for quenching (see also e.g. K ang et ani2005l: ICattaneo et alj 
i2006i:jde Lucia & Blaizot 2007j) . 

^ http://www.mpa-garching.mpg.de/NumCos/CR/Download/index.html 
^ http://galaxy-catalogue.dur.ac. uk:8080/galform/ 
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Fig. 5. — Qualitative illustration of galaxy growth and quenching in three different basic models: a "merger" model, in which systems are quenched (for any 
reason) after a major, gas-rich merger; a "halo quenching" model, in which systems are uniformly quenched when their halo reaches a critical mass Mq and 
establishes a "hot halo" gas accretion mode; and a "secular" model, in which internal galactic processes (e.g. instabilities) determine and color, independent of 
external processes. In all three models, star formation and accretion move systems to larger galaxy and halo masses in the blue cloud (blue shaded regions), and 
dry mergers move systems to larger masses in the red sequence (red shaded regions). However, the division in this galaxy-halo mass space is different in each 
case: for the "halo quenching" or "secular" cases it depends solely on halo mass or galaxy mass, respectively. In the "mergers" case, the transition line is tilted, 
as the probability of mergers depends both on galaxy and halo mass. More massive halos are more evolved, live in higher-density regions, and have more likely 
accreted other galaxies to supply a major merger, so the red fraction increases with halo mass. But at a given Mhaio, mergers are more efficient for high-mass 
systems (and initial capture more likely), so the red fraction increases with galaxy mass. Note that for all of these, we are explicitly focused on central galaxies, 
and ignore processes that may redden satellites. 
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Fig. 6. — As Figure|5] but showing the predictions from full cosmological models (again, for central galaxies only). Galaxies are color-coded by whether or not 
each model p redicts they should be in the blue cloud or red sequence. Left: Our full merger model Monte Carlo predictions. Center: The semi-analytic model of 
ICroton et alj 12006), which implements a standard halo quenching model (albeit requiring the presence of a relatively massive BH to maintain quenching). Note 
the apparent relatively low number of massive gala xies/halos owes to the sampling density of the model in its public release. Right: The modified semi-analytic 
model of Bower et al. (2006), as described in § 13.21 where we assume the strong secular (disk instability) mode that dominates the morphological transformation 
and gas exhaustion of most disks (in the model) also determines whether or not galaxies are quenched. Dashed lines in each qualitatively divide the red and 
blue populations, as in Figure[5] Despite the considerably complexity added to these models, their qualitative behavior in the Mgij— Mi,;jo plane reflects the key 
distinctions of each corresponding toy model in Figure|5] 

even without mergers, owing to the disk instability mode of 
growth). Fo r our purpose s , there fore, it is effectively equiv- 
alent to the ICroton et al.l (1200 6*) and other halo quenching 
models. But, given the strong secular mode assumed in the 
model, we easily can use it to construct an mock example of a 
semi-analytic model in which secular processes dominate the 
quenching itself. 

We do so by adopting the lBower et alj (12006 !) model, but in- 
stead of using their criterion for quenching (namely, the pres- 
ence of a hot halo), simply assume that systems which un- 



The lBower et al.l ( l2006l) models implement a strong disk in- 
stability (secular) mode, which dominates black hole growth 
and bulge formation at all redshifts, with mergers typically 
contributing only ^ 0.1% to the spheroid mass budget. How- 
ever, in the model, it is still assumed that cooling can only 
be halted in a quasi-static hot halo, and effectively galaxies 
are quenched upon crossing the appropriate halo mass thresh- 
old (like other models, the presence of a moderate-mass BH 
is technically required, but essentially all systems with suffi- 
ciently massive halos easily host a BH of the necessary mass. 
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dergo a sufficiently massive disk instability that destroys the 
entire disk will "quench." The disk stability is estimated ac- 
cording to the assumptions of the original model, based on 
e.g. disk angular momentum, scale lengths, masses, and con- 
centrations. We specifically adopt a mass threshold for the 
instability of > 2 x 1O"*M0 (i.e. assume systems in which 
less of the galaxy mass participates in the instability will not 
automatically quench, since almost all galaxies in the model 
have at least some very small mass added to the bulge via 
instabilities). We choose this value because it gives a good 
match to the total observed mass density of passive galax- 
ies and globally-averaged quenched fractions as a function 
of Mgai, but note that our comparisons are all qualitatively 
unchanged regardless of exactly how we choose the quench- 
ing criterion . We subsequently refer to this as a "modified 
[Bower et_alj (|2006|)" model, and emphasize that we are not 
plotting the predictions of the original model ( which are, for 
our p urposes, equivalent to the predictions of ICroton et al.l 
120061) . but using it to represent the predictions of a cosmo- 
logical model for secular evolution, in the case where that 
evolution dominates galaxy quenching. 

We extract the z = predictions of both models, and clas- 
sify galaxies as either red or blue following the criteria of 
the authors (namely colors (U-B) > 0.8 being "red"), al- 
though it does not change our qualitative comparisons if we 
ad opt a magnitude- dependent color limit (although, as noted 
bv lWeinmann et aLi r2006b, this reveals that high-mass galax- 
ies in both models are "too blue"). We extract these properties 
only for central galaxies - both semi-analytic models invoke 
alternative physical mechanisms such as ram pressure strip- 
ping to rapidly redden essentially all satellite galaxies. While 
there is no doubt this is an important mechanism, it has noth- 
ing to do with the models we wish to compare, and would only 
confuse the comparison we wish to highlight (and obscure the 
important differences between models). The position of these 
systems in the Mgai-Mhaio space is shown in Figure |6] in the 
same manner as Figure|5] 

To compare to these models in more detail, we construct 
a realistic Monte Carlo population of galaxies of differ- 
ent masses in different mass halos, from our merger-driven 
model. Beginning with a small halo at high redshift, hosting 
a (initially) disk-dominated galaxy (in the absence of merg- 
ers), we integrate forward in time. The average halo mass ac- 
cretion history in a ACDM universe is well-defined (here we 
adopt the av erage progenitor mass as a function of time from 
iNeistein eT al. 2006). At each point in time, the average mass 
of a disk galaxy in such a halo can be estim ated empirically , 
either from halo occupation models (e .g., .Yang et al.l 120031: 
IConrov et all 120061: IWang et al.1 l2006h . adopting the bary- 
onic Tully-Fisher relation (assuming the disk circular veloc- 
ity traces maximal halo circular velocit y, e.g . McGaug h et al.l 
12000; McGaugh 2005; B ell & de Jong| |2001), or assuming a 
constant baryon fraction in the galaxy. We henceforth adopt 
the baryonic Tully-Fisher expectation for Mdisk(Mhaio), which 
we assume does not evolve with redshift (as suggested by a 
large nu mber of observat i ons at least to z ^ 1 .5, and by some 
to z > 3; IConseHce et al. 2005; Flores et al. 2006; Bel l et all 
I2006bl: iKlssin et al . 2007; van Dokkum et al. 2004), but we 
have tried all three estimators, and find similar results. This 
is not surprising, since, as discussed in § |2] observations find 
there is little or no evolution in most general halo occupa- 
tion statistics of star-forming galaxies (i.e. average baryonic 
mass ho sted by an "un-quenched" halo of a given mass) even 
to z ~ 4 (lYan et al.ll2003[lHeymans et al.ll2006bHConroy et all 



l2006ll2007h . 

At each time, we probabilistically increase the disk 
mass with the halo mass, such that an ensemble of these 
Monte Carlo simulations always has the appropriate mean 
A^disk(A^haio) and obscrvationally measured scatter about this 
quantity. Then, we calculate the probability of a major gas- 
rich merger, specifically the probability both that the halo has 
accreted another halo hosting a galaxy of comparable mass 
(mass ratio < 3 : 1) and that the two will merge in the given 
timestep. This calculation is identical to that in § |2] (see Pa- 
per I for details), where the former probability has been de- 
termined from dark-matter simulations (i.e. the probability of 
hosting or accreting a subhalo of the appropriate mass range) 
and the latter is the ratio dt / fmereer (where fmerser is the merger 
timescale as in § |2l we generally adopt the dynamical fric- 
tion timescale in what follows, but our results are qualitatively 
similar regardless of this choice). Based on this probability, it 
is randomly determined whether or not the galaxies merge. If 
so, the final stellar mass is just the sum of the two pre-merger 
baryonic masses, and we assume zero further growth through 
star formation (although growth via dry mergers is allowed). 

We technically integrate this model only from z~10toz = 
(or where Mhaio > 10^ M©), but find the results are reasonably 
converged with respect to this choice (although in principle 
every halo may have a major merger if we integrated to infi- 
nite redshift or Mhaio = 0, these mergers are meaningless for 
our purposes as there is no significant galaxy formed inside 
the halo). Running a large sample of Monte Carlo realiza- 
tions for each Mhaio, we obtain a bivariate z = distribution of 
early and late-type galaxies in Mgai and Mhaio which reflects 
our models. The resulting predictions are shown in Figure|6] 

Although this in some sense serves as a crude toy semi- 
analytic model, we adopt this approach specifically to mini- 
mize the uncertainty owing to choices such as the modeling 
of star formation and accretion in galactic disks. Instead, we 
adopt as much as possible in a purely empirical fashion, to iso- 
late the predictions of a merger-driven quenching model (and 
not confuse these with degeneracies in modeling disk forma- 
tion). Since mergers will efficiently convert gas to stars, and 
their gravitational processes are not changed by the ratio of 
gas to stellar mass, our results are also entirely independent 
of the star formation histories in the disks - we only need to 
inform our predictions with a rough estimate of the masses 
of disks hosted by halos of a given Mhaio- Ultimately, adding 
the complications of our Mont e Carlo te sts allows us to con- 
struct a c omparison to the Cro ton et alj (2006) and modified 
iBower et al.l (1200 6) models, but yields a qualitatively simi- 
lar result to our naive cartoon expectation in Figure |5] (Note 
that there are some differences in the low-mass star-forming 
galaxies between the various models, owing to their treatment 
of star formation, but this is unimportant for any of our con- 
clusions.) 

Quantitatively, we can now integrate the results of Figure|6] 
and predict the red (i.e. merger remnant) fraction as a bivari- 
ate function of Mgai and Mhaioi Figure [T] shows this. In or- 
der to represent the real observations, we add the appropriate 
observational errors in both Mgai and Mhaio (c m,,..! ~ 0.2 dex, 
iMh,.!,, ~ 0.4 dex), fo r both our model and the [Croton et alJ 
(20M; [Bower et all ([2006[) models. This does not quahta- 
tively change any of the results, but does smooth some of the 
dependencies (and tends to remove unphysical features in the 
models caused by undersampling). 

In a global sense, the trends appear to be reasonably ac- 
curate - they agree well with the observed fraction of red 
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Fig. 7. — Top: Local fraction of red/early-type (major merger remnant) 
central galaxies as a function of halo mass, from our prediction in Figure |6] 
Linestyles adopt different estimations of tiie merger rate, as in Figure m Solid 
blue line shows the mean fraction, upper and lower green lines the fraction in 
the higher and lower stellar mass halves at each Mhaio, respectively. Bottom: 
Same, as a function of galaxy stellar mass. Green lines in this case are as the 
blue lines, but adopt a different halo occupation fit (as the dashed black fines 
in Figure[T). Black squ ares show t he ob served early-type galaxy fraction as 
a function of mass from lBell et alj 120031) . 
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Fig. 8. — Mean central galaxy baryonic mass as a function of halo mass 
(blue lines, as in Figure [T), from our prediction in Figure |6] Dotted orange 
line corresponds to the universal baryon fraction. Points sfiow the obser- 
vationally estimated m ean central galaxy stellar mass as a function of halo 
mass from iWang et aU <2006l) (HOD fitting) and Mandelbau m et aU {2003) 
(weak lensing). The merger-driven quenching model naturally predicts the 
red fraction as a function of mass and the turnover in the Mgai(Miiaio) relation 
(equivalently, turnover in galaxy M/L ratios above ~ M»), without any input 
parameters describing a preferred mass scale. 



galaxies as a function of galaxy stellar mass (iBell et alJ2003h . 
Figure [8] compares the mean Mgai at each Mhaio predicted 
from this model. Quenching associated with major mergers 
naturally predicts the turnover in Mgai(Mhaio) around Mgai ^ 
10" M0. We emphasize that there are no parameters in our 
model which have been tuned or otherwise adjusted to give 
this result - unlike halo quenching models which empirically 
adopt a specific quenching mass, we have no input param- 
eter which fixes this mass. Rather, the turnover arises self- 
consistently, as the result of major, gas-rich mergers first be- 
coming efficient at these masses, and subsequent star forma- 
tion being quenched. 

We now examine the predicted red fractions in greater de- 
tail, by breaking them down as a bivariate function of both 
Mgai and Mhaio- Figure |9] shows this, for each of the three 
models as in Figure |6] and several observational determina- 
tions. Specifically, we calculate the red/early-type fractions 
predicted as a function of Mhaio, in bins of galaxy stellar mass 
Mgai. Note that a detailed quantitative comparison with the 
observations is difficult and beyond the scope of this paper, as 
the exact absolute values of /red depend sensitively on the se- 
lection method and c onversion between g roup properties and 
halo mass (see, e.g. ICooper et alJl2006b() . But the q uaUta- 
tive trends are robust to thes e effects (see e.g.lWeinman n et al] 
2006b': ICooper e t al."2006a b). For all the model predictions 
and the observational analogues, we consider only the red 
fraction of central galaxies. In most models, satellite galaxies 
are uniformly (or close to uniformly) red, so considering the 
total (central-i-satellite) red fraction mixes the consequences 
of the physics causing quenching (what makes central galax- 
ies red) with the estimated satellite fraction as a function of 
Mgai and Mhaio (which, while importantly informing models, 
contains no information about the physics of central galaxy 
quenching). 

From the observe d group c atalogues of Weinmann et alJ 
(2006a) and Marti nez & Mur iel (2006); Martinez et al] 
(2006), which consider the same (again, for central galaxies 
only), there are a few important qualitative trends. These 
include: (1) a strong dependence of red fraction on halo 
mass, but (2) a significant residual dependence on galaxy 
mass/luminosity, (3) a lack of any sharp characteristic scale 
in Mhaio, (4) a relatively high red fraction (/red > 0.5) for 
the most massive/luminous systems even at low halo masses 
(Mhaio < lO'^M©), and (5) a similar, relatively high red 
fraction (/red ^ 0.5) for the least massive/luminous systems 
at high halo masses (Mhaio > lO'-^M©). 

In contrast to the observed trends, the ICroton et al.l (l2006h 
model is, as expected, similar to a pure halo quenching model 
- there is a sharp transition from uniformly low red fractions 
(/red ^0.1) below the halo quenching mass (-a few IO'^Mq) 
to uniformly high red fractions above this halo mass, with a 
weak residual dependence on galaxy mass. The low red frac- 
tion at small halo masses also forces these models to assume 
a high red satellite fraction at these masses (in order to match 
the global red galaxy mass functions), in disagreement with 
observations (Weinmann et al. 2006^. 

The public (original) version of the iBower et alJ (l2006h 
model yields an essentially identical prediction to the 
ICroton et al.l (l2006l) model in this space, as the development 
of a hot halo is assumed to be the key criterion for quenching. 
The modified Bower et al. (2006) model which we consider, 
on the other hand, is quite similar to a pure secular model 
(as expected), with /ed nearly independent of Mhaio at each 
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Fig. 9. — Red/early-type fraction /^d (of central galaxies only) as a bivariate function of stellar mass/luminosity and halo mass/local environment. We 
specifically exclude satellites, as they tend to be uniformly red (making the predicted red fractions degenerate between central galaxy quenching mechanisms and 
the satellite fraction as a function of Mhaio and Mg^i). Top Left: Observed /led of centra l galaxies as a fu nction of host halo mass (esti mated from matching grou p 
catalogues to halo ma ss functions) in bins of galaxy r-t)and magnitude, from Weinm ann et all <2006aL solid line, filled points) and lMartmez & Muriell I200(il) : 
IMartmez et all )200fl dashed line, open points). (Note there appears to be some small fraction of massive galaxies in small halos in each panel: this owes to 
scatter in the halo and stellar mass estimators, but has no effect on the conclusions.) Top Right: Predicted /i-ej of central galaxies from our merger model, as 
a function of halo mass in bins of galaxy stellar mass, as la beled (bins of a giv en color/style roughly correspond to the observed r-band absolute magnitude 
ranges of the same color/style). Bottom Left: Same, from the Croton et alj flQOQ) halo quenching model. We qualitatively label the quenching halo mass, which 
separates uniformly low and uniformly high in this model. Bottom Right: Same, from the modified Bower et al. (2006) secular model. Dotted lines show the 
Mijaio -independent red fraction for each M^^\ if the model were strictly dependent on only M^^\ (the model treats satellites and central galaxies differently, so the 
normalizations of these /.ed estimates does not agree at low Mg^i). The behavior of the three models is qualitatively different, as in Figure[6] with a merger model 
predicting a joint dependence on Mgai and Mhaio distinct from the halo quenching or secular models. 

Mgal. There is some weak dependence, because galaxies liv- 
ing in high-mass halos tend to have earlier formation times, 
meaning that their progenitor disks were more compact and 
therefore (according to the model assumptions) more prone 
to massive instabilities, but the primary dependence of /red is 
clearly on galaxy mass. 

Neither of these predictions agrees qualitatively with the 
observations. The prediction from our merger model, how- 
ever, matches these features - t he dependenc e on Mhaio is 
stronger than that in the modified iBower et alj (120061) model 
(or a "toy" secular model) in considerably better agreement 
with the observation s, but the residual de pendence on Mgai is 
stronger than that in ICroton et al.l (l2006i) . There is no sharp 
transition at some specific Mhaio, and the red fraction of mas- 
sive s ystems remain s relatively high at lower Mhaio, in contrast 
to the Croton et al.l (12006) predictions. However, there is still 
a significant dependence on halo mass, and low stellar mass 
systems do bec ome red a t large Mhaio, in contrast to the secu- 
lar/modified BowerelaD (2006) model predictions. 

Observationally, when red fractions are quantified as a 
function of quantities such as galaxy density p or surface den- 
sity E, there is some ambiguity in what these quantities rep- 
resent. To lowest order, they serve as tracers of halo mass 
and are directly co mparable to p redictions such as those in 
Figure m However. iBaldrv et al.l ( l2006l) and others have sug- 
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Fig. 10. — Dependence of red fraction on density at small scales (left) 
and large scales (right), at fixed halo mass (i.e. considering /icd/(/icd) versus 
density (1 + 5r)/ ((1 +5r)) at fixed Mhaio)- Points show the observations from 
Blanton et al. (2006), for SDSS groups with different total group luminosi- 
ties (as labeled; this should be a good proxy for total group halo mass). Lines 
show our prediction, which has an increasing /caiiy with overdensity on small 
scales, as a local galaxy overdensity implies an increased probabihty of ma- 
jor mergers. Solid fine is for z = 0, Mhaio = M» 1 .5 X lO'^ /r' Mq halos, 
dashed lines show how this changes for more massive (Mhaio ~ lO" ' ^0. 
shallower dependence on density) and less massive (Mhaio ~ lO'^/r'M©, 
steeper dependence on density) halos. The merger-driven quenching hypoth- 
esis naturally explains a dependence on small-scale overdensity, similar to 
that observed. 
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gested that the trends in /red with these quantities argue for 
some level of dependence on environment, even after account- 
ing for the primary depe ndence o n halo mass. 

In greater detail, Blant on et alj (|2006) investigate this pos- 
sibility by determining /red in SDSS groups as a function of 
local density (1 + 6/, defined as the galaxy number density 
within a radius r relative to the mean number density in that 
radius) on various scales, in narrow bins of total group lumi- 
nos ity (which should be a good proxy for gr oup halo mass, 
see lYang et ani2005l: Ivan den Bosch et al.l l2006). At fixed to- 
tal group luminosity (roughly equivalent to fixed group parent 
halo mass), they find no evidence for an additional depen- 
dence of red fraction on large-scale environment, measured at 
projected radii 6 < r < 10 and 0.3 < r < 1 /;"' Mpc. However, 
at small radii 0.1 < r < 0.3/i~' Mpc they find a significant de- 
pendence of /red on density for all group luminosities (halo 
masses) which th ey consider. A simi l ar res ult is found by 
iPark et al.1 (l2006l) and lKauffmann et al.1 (|2004 . 

In a halo quenching or secular model, this is difficult to 
explain, as quenching depends only on either halo mass or 
internal galaxy properties, respectively. However, as we have 
discussed for both ongoing/recent galaxy mergers and quasars 
in Paper 1, mergers are more likely to occur in regions with 
galaxy overdensities on very small scales. The bias to increas- 
ing fractions of merger remnants with increasing small scale 
density in Paper 1 (see Figures 7 & 17 in that paper) directly 
translates to a prediction for the dependence of /ed on small 
scale overdensity, which we show in Figure [TO] The merger 
hypothesis provides a natural explanation for the observed de- 
pendence on small-scale overdensities. 

We caution, however, that this explanation is not unique. 
If satellites are preferentially red, then a simple autocorrela- 
tion function or dependen ce of overall red frac tion on density 
(such as that observed by iBlanton et al.ll2006h of red galax- 
ies will see a similar effect (with the excess on small scales 
reflecting the abundance of satellite galaxies). Furthermore, 
mergers (by definition) consume some of the galaxies that 
(initially) define the small-scale overdensity, so it is not clear 
how much of this effect might be wiped out by the merg- 
ers themselves. In other words, seeing this effect weakly 
does not necessarily argue against a merger-driven model for 
quenching, and theoretical study in cosmological simulations 
is needed for more detailed predictions. 

A more rigorous test of this would be to compare the cross- 
correlation of central red galaxies with all other galaxies, 
i.e. to more directly test whether central red galaxies pref- 
erentially live in regions of small-scale galaxy overdensity. 
Early a nalysis along these l ines does suggest a similar con- 
clusion dMasiedi et al .1120071) . in agreement with these predic- 
tions (and difficult to reconcile in models where quenching is 
a pure function of galaxy or halo mass). 

3.3. Redshift Evolution of Quenched Fractions and 
Color-Morphology-Density Relations 

We next examine the redshift evolution of the trends in red 
galaxy fraction with stellar and/or halo mass. First, we return 
to our prediction for the mass density of passive systems as 
a function of redshift (Figure |3]l. As noted there, a merger 
quenching model predicts a (relatively high) density of pas- 
sively evolving, quenched systems in good agreement with 
that observed at high redshifts z > 2-3. This is in strong 
contrast to many pure hot-halo quenching models, in which 
cold accretion within a hot halo persists at high redshifts 
Z > 2. Figure [TTI contrasts the merger-driven prediction with 
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Fig. 1 1. — Integrated stellar mass density of red or early-type galaxies 
as a function of redshift. The predictions from a merger quenching model 
(with shaded range reflecting the range of predictions from our different 
adopted models) and observations are as in Figure |3] We compare the pre- 
dictions of a pure sec ular model (based on the observed total mass function 
in Font ana et al.l 12006). where we assume the red fraction as a function of 
galaxy mass is identical to that at z = 0), and various halo quenching models 
(where we assume all systems above the quoted halo quenching masses in 
our halo occupation model are quenched). In order to form massive galaxies 
at high redshift, the traditional halo quenching models must allow cold ac- 
cretion flows in all but the most massive halos at z > 2 - 3, yielding almost 
no red galaxies at these redshifts. In contrast, by allowing reddening to occur 
in a range of halo masses (which might otherwise continue accreting), the 
merger and secular model produce a sufficient density of passive galaxies at 
high redshift. 

that from several halo quenching models (Dekel & Birnboiml 
2006; Cattaneo et al. 2006; Croton et al. 2006), where we as- 
sume (since we are ultimately just making a qualitative com- 
parison) that all systems above their quoted quenching halo 
mass thresholds are red/pas sive (but our explicit compari- 
son with these models in § 13.21 suggests this is actually a 
good approxim ation to their predictions). Note that for the 
ICattaneo et al.l 12006), this is a comparison with the halo 
quenching-only version of their model (since they also con- 
sider a model in which, like ours herein, major merger rem- 
nants are automatically quenched). We then use our halo oc- 
cupation model to determine the "red" mass density. 

The result is similar in each halo quenching model: above 
z ~ 2-3, the densit y of passive galaxies plum mets. In de- 
tail, for example, the lDekel & Birnboiml (l2006l) derivation of 
the hot halo quenching mass predicts that at z ~ 3.5, cold 
flows continue within all Mhaio ^ 10 ''^Mq halos, which al- 
lows only a completely negligible maximum red galaxy mass 
density (since such halos are extremely rar e at high redshifts). 
As demonstrated bv lCattaneo et alj (l2006h . introducing these 
"cold flows in hot halos" is necessary for these models to 
match the overall density of massive galaxies and cosmic star 
formation rate density at high redshift. This owes to the steep 
step-function transition from unquenched to quenched sys- 
tems around the quenching mass in such models (see Figure|9|l 
- if the quenching mass is lowered (to make for more red 
galaxies), then the models will quench systems too early, and 
not form any high-redshift massive galaxies in the first place. 

In contrast, a merger driven model is able to predict that 
the appropriate fraction of these massive, high redshift galax- 
ies are passive. This is because it allows for reddening to 
be somewhat uncoupled from halo mass; i.e. systems in mas- 
sive halos might generally continue to accrete, but some frac- 
tion can redden and build up sufficient mass density of passive 
galaxies (without reddening all systems of these masses and 
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destroying the ability to make massive galaxies in short cos- 
mic times). 

In addition, we compare the simplest secular model (where 
we just adopt our halo occupation model and assume the red 
fraction as a fun ction of galaxy mass is identical to that at 
z = 0; as in § 13.2b . and (for the same reasons) find that it is also 
able to reproduce the observations. This should not be surpris- 
ing, since in this toy model, matching the overall galaxy mass 
density is implicit, and the red fractions are as high as they are 
at z = 0. However, we will show that there are other aspects 
in which the redshift evolution of red fractions predicted by a 
secular model are in conflict with the observations. 

Interestingly, simulations suggest that hot halos often de- 
velop at lower masses th an the halo qu enching models re- 
quire (- 10"-** Mq; see Kere s et aLll200 5) - in other words, 
the possibility that such halos are necessary for quenching 
is viable, since this mass threshold allows for the possibil- 
ity of sufficient populations of passive, high-redshift galaxies. 
However, the idea that such halos are sufficient for quench- 
ing (as is effectively true in the halo quenching models) is not 
viable, since it prevents the formation/growth of galaxies be- 
yond this mass threshold and cannot form sufficient numbers 
of massive galaxies nor yield a sufficiently high global star 
formation rate at high redshift. 

Figure [12] shows the mean predicted early-type fraction in 
our merger-driven model as a function of halo mass (as in 
Figure|7]i at different redshifts. We also plot this as a function 
of the estimated clustering amplitude for each Mhaio (at each 
redshift), and the (approximate) corresponding comoving cor- 
relation length ro (where we define ro = S/z"' Mpc [b(z)]^ , 
where 7 « 1.8 as measured locally). At high redshifts, sys- 
tems in only the most massive halos (i.e. most extreme over- 
densities) have sufficiently rapid merger rates that they will 
have built up large red fractions. In terms of Mhaio, the evolu- 
tion appears relatively weak, but in terms of the clustering am- 
plitude, it is more obvious (this owes to massive halos being 
rarer and corresponding to significantly higher-density peaks 
at high redshifts). 

Fig ure[T3]compares thes e predictions to recent observations 
from Cooper et alJ (l2006al ). The same trends are evident for 
both our predictions and the observations. Looking at typ- 
ical halos, expected in average regions of the universe, the 
color-density relation (more specifically, the dependence of 
/early on halo mass) nearl y vanishes by z ^ 1 .5. Similar trends 
have been observed by Gerke et alJ (12007) and Nuijten et alJ 
(j2005). This does not mean that there is no such trend - it 
does not become prominent until higher-density, rare high- 
halo mass peaks are probed, as in Figure [121 This is also 
observed - for the most massive, red galaxies at z ~ 2 - 4, 
a color-density relation is seen for dense environments with 
f() ^ 10/;~'Mpc, similar to our predictions. Again, we em- 
phasize that a detailed quantitative comparison is outside the 
scope of this paper, as it is sensitive to the galaxy and color 
selecti on method and the e xact definition of galaxy environ- 
ments (ICooper et al.l l2006a). but the qualitative trends should 
be robust. We also caution that the observations (presently) 
mix central and satellite systems, although above moderate 
~ luminosities, satellite systems should only be a small 
fraction of the observed populations. 

In Figure [14] we compare the evolution in the red fraction 
predicted as a function of halo mass and redshift from our 
merger-driven model to that predicted in a secular or halo 
quenching model. For the secular model, we assume that the 
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Fig. 12. — Top: Evolution in the predicted early-type fraction as a function 
of halo mass (see FigureQ with redshift (here for just our standard model). 
Bottom: Same, but as a function of the halo clustering amplitude (calcu- 
lated for each M^aio: z) or (rough) equivalent congelation length rg. The fig- 
ure should be interpreted as reflecting the scales on which a color-density or 
morphology-density relation will manifest at different redshifts. At high red- 
shifts, the highest-density regions (e.g. > 2, ro > 10 at z 2) will have built 
up a color-density relation, while at low redshifts the color-density relation 
will have built up even in field populations (b ~ 1). 

red fraction is a pure function of galaxy mass (calibrated at 
z = as in Figure O - so the difference in red fraction as a 
function of halo mass reflects only evolution in the typical 
stellar masses hosted in halos of a given mass. As we have 
noted (see §|2]i, this evolves relatively weakly, and as a conse- 
quence there is little evolution in the red fraction as a function 
of halo mass in secular models. It is clear in Figure [T4l that 
(in the secular model) the color-density relation for "typical" 
^ 10" - 10'^/i~'Mq mass halos does not vanish at z ^ 1.5 
as is observed. Although there might be some apparent evo- 
lution (since, as noted in Figure [121 the clustering of halos 
of fixed mass will vary with redshift), there is little evolu- 
tion in the color-density relations in terms of halo mass in this 
model, and it is not obvious how any fundamentally secular- 
dominated mod el can avoid this pre diction. 

Interestingly, Cassata et alJ (l2007h observationally estimate 
that the pseudobulge fraction as a function of redshift may 
exhibit this behavior (being essentially constant for a given 
galaxy/halo mass), as expected if these are formed via secu- 
lar mechanisms, and ISheth et aP (|2003) find similar results 
for th e evolution in disk bar fractions. But Ca ssata et alj 
( l2007h (and others) estimate that this pseudobulge popula- 
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Fig. 13. — As Figure nn but comparing with observations. Top Left: Pre- 
dicted red fraction versus lialo mass and redsliift, for typical lialo masses (neai' 
M* ) at z ~ — 1 . Bottom Left: Observed red fraction verses densi ty (which 
should trace A/haio. to lowest order) at different redshifts from Cooper et al.l 
(20063). The predicted and observed relations flatten with redshift in a sim- 
ilar manner. Top Right: Predicted red fraction versus redshift, for different 
halo masses (log(Mi,aio//i"' Mq) = 10.0, etc., as labeled). The most massive 
halos are in high-density regions which evolve most rapidly at high r edshifts. 
Bottom Right: The observed red fraction versus redshift from Co oper et al.l 
(2006a) for the highest-density and lowest-density 1/3 of systems at each 
redshift (shaded ranges). An exact quantitative comparison is sensitive to 
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Fig. 14. — R e d or quenched fraction as afunction of halo mass andred.shift 
(as Figures rm&fm. for a pure secular model (top) and a pure halo quenching 
model (bottom). Because the average total mass of a disk in a halo of a 
given mass evolves weakly with redshift, there is little evolution in the color- 
density (morphology-density) relations in the secular model, in contrast to 
the observations. 



tion accounts for only < 5% of massive spheroids. Future 
observations and d irect comparison with e.g. the model of 
iBower et all (l2006h can place stronger constraints on these 
distinctions, but even at present, the observations sug gesting 
evolution in the color-density relations at z > 1 dCooper et alJ 
|2b06a; Gerke et al. 2007; Nuijten et al. 2005') appear to con- 
tradict the basic prediction of a model in which secular pro- 
cesses dominate red galaxy formation. 

We also consider the predictions from a halo quenching 
model in Figure [14] Specificall y, we adopt the red fraction 
as a function of halo mass from lCroton et al.l (l2006l) at z = 
(a near step-function rise near the halo quenching mass), and 
renormalize it by the evolution in the halo quenching mass 
with redshift (i.e. shift the step function to whatever halo mass 
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Fig. 15. — Average predicted major merger rate for central galaxies in 
halos of a given z = mass (as labeled). Lines show the average rate of 
disk-disk (gas-rich), mixed morphology (disk-spheroid, but also gas-rich), 
and spheroid-spheroid (gas-poor, i.e. dissipationless or "dry") mergers. In 
an integrated sense, mixed-morphology mergers are a relatively small con- 
tribution to global merger rates. Dissipationless mergers, however, become 
dominant in the massive (Mi,aio > lO" Mq) systems at late times (once most 
systems of the given mass have already undergone at least one major merger). 

corresponds to the model halo quenching mass at that red- 
shift). The predictions are somewhat different from those of 
our merger-driven model, but in this case that owes mostly 
to the difference in red fraction as a function of halo mass 
at fixed redshift (as in Figure |9]l. Future observations of the 
red fraction as a bivariate function of galaxy and halo mass 
at different redshifts can break these degeneracies, but for 
now we note that the evolution is also qualitatively consistent 
with the observations - in both this model and the merger- 
driven model, the field color-density relation begins to disap- 
pear around z > 1 , with the buildup of the color-density rela- 
tion occurring at higher masses (higher density environments) 
at higher redshifts. 

3.4. The Role of Dissipationless or "Dry" Mergers 

Our model directly yields the major merger history of a 
given galaxy population. We therefore briefly quantify this 
as a function of galaxy properties, decomposing the types of 
mergers in comparison to various observations. Figure [15] 
for example, shows the average major merger rate predicted 
by our model (specifically our Monte Carlo realization de- 
scribed in § 13.2b for central galaxies of a given mass, decom- 
posed into the average rates of disk-disk mergers (i.e. two sys- 
tems which have not yet undergone a major merger), mixed- 
morphology mergers, and spheroid-spheroid mergers (i.e. two 
systems which have both undergone previous major mergers). 

Once the fraction of systems which have undergone at least 
one major merger at a given Mg^i and Mhaio becomes large, 
spheroid-spheroid mergers will naturally become the domi- 
nant type of merger. In other words, merger efficiencies are 
not especially sensitive to galaxy types (at fixed mass), and so 
reflect the abundance of merged or un-merged systems. We 
note that only this category, the spheroid-spheroid mergers, 
will be (in our model) gas-poor, dissipationless or "dry" merg- 
ers. We show the results just from our default model, but note 
that they are qualitatively similar regardless of our choices of 
halo occupation models, subhalo mass functions, or merger 
timescales. 

The rate of dry mergers in massive systems is consistent 
with observational estimates dBell et al.ll2006al : Ivan DokkumI 
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Fig. 16. — Top: Average number (with ~ 1 cr range shown as error bars) 
of major mergers in the history of a z = early-type galaxy (i.e. galaxy with 
at least one major merger in the past) as a function of galaxy stellar mass 
(left) or host halo mass (right). Bottom: Fraction of early-type galaxies (as 
a function of galaxy stellar mass [left] or host halo mass [right]) which have 
undergone just their one, initial (gas-rich) spheroid-forming major merger 
(blue), or more than one major merger (red). Orange lines decompose the 
red line into systems which have undergone exactly 2 major mergers (i.e. 1 
gas-rich, and generally 1 spheroid-spheroid or "dry" merger; dot-dashed), 3 
major mergers (2 dry mergers; dashed), and > 4 major mergers (> 3 dry 
mergers; dotted). 

l2005h of roughly ^ 0.5- 1 dry mergers per massive elliptical 
smce z ~ 1 (i.e. - O.lGyr-i in Mhaio > lO'^M© halos at z < 
1). Although briefly important in the transition between the 
dominance of gas-rich and gas-poor mergers, mixed morphol- 
ogy mergers are an intermediate phenomenon - most galaxies 
that have undergone only their initial, gas-rich, spheroid form- 
ing major merger were produced in disk-disk mergers, and the 
evolution of massive systems with multiple mergers is domi- 
nated (at late times) by spheroid-spheroid mergers. Note that 
for most of our predictions above, only the fact that the merger 
remnant is quenched is important, although the morphologies 
of the progenitors can change the "type" of merger. 

Integrating these rates to z = 0, Figure [16] shows the mean 
number of major mergers (and fraction of spheroids with a 
given number of previous major mergers) as a function of 
spheroid or host halo mass. We show this only for spheroids, 
since (by definition) the fraction of late-type systems (i.e. sys- 
tems which have not undergone a major merger) is identical to 
our blue fractions in § 13.21 There is a general trend for more 
massive systems to have experienced a larger number of ma- 
jor mergers, as expected since such systems form earlier and 
in more dense environments. The trend is somewhat steeper 
as a function of galaxy mass than as a function of halo mass 
- this is also expected, since at a given Mhaio, a larger number 
of mergers builds a more massive galaxy, steepening the trend 
in number of mergers as a function of Mgai. 

The number of mergers as a function of mass is simi- 
lar to that predicted by vario us semi-analytic models (e.g., 
Khochfar & Burkert 2003; Cr oton et all 120061: iKang et alj 
20051) and relatable to (although not identical to, owing to the 
differen ce in definitions) the "effective number of progeni- 
tors" in iDe Lucia et al.l (l2006l) . This is expected, as mergers 
are a dynamical inevitability. There might be some differ- 
ences owing to various prescriptions for merger timescales 
or different populations of galaxies in a given halo, but the 
general results of Figure [16] are robust. The nature of these 
mergers, however (whether they are, for example, gas-rich or 



gas-poor), does depend on the model. 

It is well-established that there is a general dichotomy in the 
properties o f elliptical galaxies (e.g . [Bender et al.'1992','1993l; 
Kormendy 1977t iKormendv etalj r2007; Laueretal. 2006^ 
Whether or not the division is strict (see e.g. iFerrarese et aP 
.2006) . the most massive ellipticals tend to be slowly rotating, 
anisotropic systems with boxy isophotal shapes and central 
core profile deviations from a pure Sersic profile. Less mas- 
sive ellipticals, including most of the ^ L^, population, tend to 
be more rapidly rotating, with disky isophotal shapes and cen- 
tral light cusps. The transition between the two occurs at ap- 
proximately My ^ -21 .5, and is commonly thought to derive 
from the difference between systems which have undergone 
just their initial, spheroid forming and gas-rich merger (which 
will dominate the less massive systems) and those which have 
undergone subsequent spheroid-spheroid dry mergers (which 
will dominate the most massive systems). 

Indeed, detailed numerical simulations have shown that 
gas-rich disk mergers, and only gas-rich mergers, repro- 
duce the detailed distributions of kinematic properties of 
the less massiv e/rapidly rotating /disky/cuspy ^ L^, ellipti- 
cal population iCox et all l2006bl) - including their rotation 
properties, kinematic misalignments, i sophotal shapes, el - 
lipticities, central light cusps (Mihos & Hernquisj Il994ah . 
velocity profiles (Naabetal. 2006a), kinematic subsys- 
t ems (Hernqui st & Barnes 1991), and internal correlations 
([Robertson et 31.112006 ). Likewise, mergers of kinematically 
hot systems, i.e. spheroid-spheroid mergers, with little gas 
content, are required to produce the combination of boxy 
isophotal shapes, anisotropy, and low rotation seen in most 
massive ellipticals (Naabetal. 2006b), and the commonly 
adopted theory of "scouring" by a binary black hole in the 
formation of central cores also requires that the mergers have 
very little cold gas content (since even ^ 1% of the stellar 
mass in cold gas falling to the center is 3> Mbh, and would 
allow for rapid coalescence of a merging binary). Mergers of 
disks, even when gas-free, cannot reproduce the combination 
of low ellip ticities and little rotation seen in the most massive 
spheroids ( Cox et al.ll2006bl) . 

Figure [T7[ therefore compares the fraction of cusp/core, 
disky/boxy, and rotating/isotropic ellipticals as a function of 
galaxy stellar mass to our estimate of the fraction of z = sys- 
tems for which the last major merger was a gas-rich, spheroid- 
forming merger, or for which the last merger was a (subse- 
quent) spheroid-spheroid dry merger The agreement is good, 
for all three indicators. Both the trend in the fraction as a func- 
tion of mass, and the transition at Mgai ^ 2-3 x 10" M0 are 
predicted by our model. This transition point is robust, with 
a rough ~ 0.2 dex systematic uncertainty owing to the exact 
version of our model which is used to calculate the merger 
histories (within the range of uncertainties from the observa- 
tions). 

This additionally puts strong constraints on other models, 
shown in Figure [18] In a pure secular model, most ellipti- 
cals are formed via disk instabilities - this is already in con- 
flict with the kinematic arguments in §[T]( which find that disk 
instabilities generically form pseudobulges, not the classical 
bulges that dominate the spheroid population at the masses 
of interest here), but in addition, these systems therefore are 
already gas exhausted and quenched by the time they un- 
dergo their first major merger Nearly all major mergers in 
such a model, then, constitute dry, spheroid-spheroid merg- 
ers. Adopting our calculated merger histories as a function 
of mass (which are not sensitive to our model for the colors 
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Fig. 17. — Predicted fraction of splieroids for wliicli tlie last merger 
was a gas-rich, spheroid-forming major merger (blue), or for which the last 
merger was a subsequent spheroid- spheroid (dry) major merger (red), as a 
function of galaxy stellar mass. We compare with the observed fraction of 
cusped or cored-central profile ellipticals ( Lauer et al. 2006), the fraction 
of ell i pticals with disky (04 /o > 0) or boxy isophotal shapes IBender et alj 
119931 : IPasquali et al J 120071 circles and stars, respectively), and the frac- 
tion of rapidly (log(v/a)* > -0.15) or slowly rotating/isotropic ellipticals 
(Sender et al. 1992). The dichotomy between elliptical types is reproduced 
well, if dry mergers form cored, boxy, slowly rotating remnants (as suggested 
by numerical simulations). In each panel, the solid fines are tfie predictions 
of our merger model, dotted and dashed lines show th e pr edictions of secular 
and halo quenching models, respectively (see Figure [l8l for clarity just the 
dry merger fractions are shown). 

and morphologies of the merging systems), but using the pure 
secular model criteria (as in Figure |9]l to determine whether 
the progenitor galaxies are already red (i.e. whether or not 
the merger is dry), we obtain the prediction that the transition 
to dominance of dry, spheroid-spheroid mergers should occur 
at masses an order of magnitude lower, Mgai ~ 10'"Mq, an 
order-of-magnitude contradiction with the observations. 




Fig. 18. — Predicted fraction of spheroids for which the last merger was a 
gas-rich, spheroid-forming major merger (blue), or for which the last merger 
was a subsequent spheroid-spheroi d (d ry) major merger (red), as a function 
of galaxy stellar mass as in Figure [TTl but for different models. The predic- 
tions of our merger model (solid) are compared to those from a pure secular 
(dotted) or pure halo quenching (dashed) models. Secular models quench 
and form bulges via disk instability, so most mergers even at low mass are 
dry/spheroid-spheroid; halo quenching models do not prevent spheroids from 
re-forming disks below the halo quenching mass, so only the most massive 
mergers are dry/spheroid-spheroid. Both predict a transition point between 
cuspy/disky/rapidly rotating and cored/boxy/slow r otat ing ellipticals an order 
of magnitude discrepant from that shown in Figure flTl 

Likewise, the simplest pure halo quenching models fail to 
reproduce the observed elliptical dichotomy. In such models, 
a substantial fraction of galaxies will experience their first ma- 
jor merger before the system crosses the quenching halo mass 
threshold. Since they are below this threshold, the system 
will re-accrete gas and re-form a disk, even for major merg- 
ers occurring just a short time (~ a couple Gyr) before the 
halo grows sufficiently massive to cross this threshold. The 
next major merger will therefore be assumed (in the model) to 
be a gas-rich disk merger, instead of a dry spheroid- spheroid 
merger. 

We again calculate the effects of this in Figure [18] using 
our merger histories but assuming that all syst ems b elow the 
halo quenching mass threshold of Croton et al. (l2006l) (at each 
redshift) will re-accrete and remain gas-rich. The result is 
that only the most extremely massive systems, which crossed 
the halo quenching mass threshold at early times, have had 
sufficient time to then undergo subsequent multiple mergers 
(yielding at least one spheroid-spheroid dry merger). We ob- 
tain the prediction that the transition point between elliptical 
types should occur at masses an order of magnitude higher 
than in our merger-driven model, at Mgai > IO'^Mq, once 
again an order-of-magnitude contradiction with the observa- 
tions. 

This is further demonstrated in the recent analysis by 
iKang et al] (l2007l) . who consider the predictions for the num- 
ber of boxy ellipticals in a similar halo quenching model. 
For the reasons given above, their model predicts that a neg- 
ligible (<C 10%) fraction of early-type galaxies have under- 
gone a true dry (spheroid-spheroid) merger, at all masses 
Mgai < 1O^^M0. In order to match the observations in Fig- 
ure [TtI they are forced to assume that any major merger 
with a gas fraction /gas < 0.1 produces a boxy elliptical. As 
a consequence, such a model predicts that ^ 1/3-1/2 of 
"boxy" systems are actually formed in mergers of two disk 
galaxies (low gas-fraction. Milky Way-like disks), with the 
remaining ^ 1/2-2/3 formed in what we would call gas- 
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rich, mixed morphology mergers. However, numerical sim- 
ulations of major mergers of disk-dominated galaxies with 
gas fractions /gas ^0.1, an d kinematic a nalysis of compara- 
ble local merger remnants dRothberg & Joseph 2006a), have 
clearly established that such mergers do not, in fact, generi- 
cally produce boxy elliptica ls. They instead pr oduce systems 
resembling disky ellipticals (iNaab et al.l2006al) . with substan- 
tial central c usps (th e cusps do not disappear even at low /gas; 
Mihos & Hernquist 1994ah . and high ellipticities jCox et al.l 
2006bl) . 

If we instead begin knowing the properties of mergers that 
form different types of ellipticals, the observations lead us to 
conclude that some form of quenching must be able to oper- 
ate, at least temporarily, in massive spheroids after their first 
formation epoch, in order that they be truly dry/spheroidal 
in subsequ ent mergers. A similar conclusion is reached by 
iNaab et all (l2006bl) . who find, using numerical simulations to 
measure the distribution of spheroid isotropy/anisotropy that 
would be observed in a remnant of a merger of specific types 
of progenitor galaxies, that matching the trend and transition 
to the dominance of anisotropic galaxies requires the quench- 
ing of all systems with massive bulges (> 3 x lO'^M©, in their 
case), effectively identical to our merger quenching criterion. 
This demonstrates the strong constraints that can be placed 
on models for how systems quench and become red galaxies, 
given the specific kinds of galaxy mergers required to produce 
the correct distribution of elliptical kinematic properties as a 
function of mass. A more detailed investigation combining 
these cosmological predictions with detailed numerical simu- 
lations, to study the effect of such mergers on the kinematics, 
internal correlations, and redshift evolution of massive ellip- 
ticals is outside the scope of this paper, but is an important 
subject of future work. 

4. THE PHYSICS OF QUENCHING 

We now turn to a discussion of the physical mechanisms by 
which mergers might both terminate significant star forma- 
tion, and result in a system which can maintain relatively low 
star formation rates. It is not our intent in this discussion to 
prove that a particular mode of feedback, for example, must 
quench subsequent star formation, but rather to highlight the 
physical processes that operate in mergers and their possible 
or likely effects on the intergalactic medium (IGM) and sub- 
sequent cooling of halo gas. 

To do so, we examine a large suite of hydrodynamical sim- 
ulations of disk galaxies and major mergers between them, 
described in detail in iRobertson et al.l ( 12006 ). The simula- 
tions are high-resolution (spatial resolution ~ 20 pc in the best 
cases), fully hydrodynamic calculations which incorporate a 
self-consistent, observationally motivated model for a mul- 
tiphase interstellar medium, star formation, supernova feed- 
back, and black hole accret i on and feedback (for details, see 
ISpringel & HernquistillOOl ISprineel et al.]l2005bl) . We con- 
struct stable, equilibrium disk gala xies to either merge o r 
evolve in isolation as described in IRobertson et al.l (l2006 l). 
with a dark matter halo, gas and stellar disk, and bulge compo- 
nent relevant for observed galaxies of the given mass and red- 
shift, with e.g. the scale length of these components set by the 
appropriate concentration and spin parameter as a function of 
mass and redshift. We specifically consider a subset of disks 
with baryonic masses Mgai « 10'°, 10", IO'^Mq, and initial 
simulation gas fractions of /gas = 0.4 and 0.2. In our merger 
simulations, we place two identical disks with a relative incli- 
nation of ^ 60° (representative of most random encounters) 



on a parabolic orbit and allow the system to evolve until it has 
completely relaxed (usually ^ 2.5-3Gyr after the merger). 
The simulations were performed usin g the code Gadget-2 
(Springel 2005), a fully conservative dSpringel & HernquistI 
120021) implementation of smoothed particle hydrodynamics 
(SPH). 

We show just the results from these cases in what follows, 
for simplicity, but n ote that we have s urveye d a much wider 
parame t er space in IRobertson et al.l (l2006l) : iHopkins et al.l 
( l2006ah : ICox et all (12006b'), varving masses fromMgai 10**- 
1O'^M0, gas fractions /gas = 0.05- 1, concentrations, bulge- 
to-disk ratios, and (in mergers) orbital parameters, relative 
disk inclinations, and merger mass ratios. We ultimately find 
qualitatively similar results in all these cases, and for our pur- 
poses the subset of simulations shown is representative of the 
important qualitative effects. 

In each simulation, we assume any stars present at the be- 
ginning were formed according to the best-fit observed t- 
m odel star format i on his tory for disks of the given mass 
in iBell & de Jond ( |2000|) (appropriate for the redshift at 
which the simulation is initialized). The stars formed dur- 
ing the simulation have ages and metallicities determined self- 
consistently. Knowing the star formation and enrichment his- 
tory of all stars in the simulation, we integrate to calculate 
the mean (B - V) color of the galaxies at each time, using 
the stellar population synthesis mod els of iBruzual & Charlon 
(l2003h with an assumed IChabrieii (l2003h IMF (similar in 
the predicted colors to our generally adopted diet Salpeter 
IMF). Because the simulation also includes gas, we can self- 
consistently integrate along the line of sight to all star parti- 
cles and calcul ate the appropr i ate du st reddening and extinc- 
tion, following iHopkins et alJ (l2005h . However, because we 
are primarily interested in times after the merger (i.e. after 
most gas is exhausted in star formation) and average trends, 
we find this makes little difference. 

4.1. "Transition" : Termination of Star Formation in Major 

Mergers 

It is relatively easy to see how a major merger can terminate 
star formation in an immediate sense. The rapid consumption 
of gas in the final stages of the merger, potentially coupled 
with expulsion by feedback mechanisms, allows for a sharp 
truncation in star formation. Figure[T9]illustrates this. 

We first consider simulations of "truncated" disks; i.e. disks 
which are completely cut off from a gas accretion supply. We 
construct appropriate disks of the masses and gas fractions 
shown in Figure[T9]and evolve them in isolation (allowing no 
further gas accretion). Technically, in terms of the stellar pop- 
ulations and disk properties, the plot assumes that disk accre- 
tion is truncated at z = 2, appropriate for most of the star for- 
mation in present-day early type galaxies (e.g.'Gallaz zi et al] 
2006, and references therein), but the qualitative result is al- 
most identical regardless of when we initialize the simulation. 
The star formation rate, plotted as a fraction of that at the 
onset of the simulation (since the optical colors here are pri- 
marily influenced by the relative decline in star formation), 
decays weakly. In fact, this drop is similar to that expected in 
the simplest models. For any disk which obeys a r-model star 
formation history M cx exp(-f/r) prior to truncation, and a 
Kennicutt-Schmidt star formation law Esf oc Sga^ (iKennicuttI 
Il998h . it is straightforward to calculate the subsequent evo- 
lution in the star formation rate if the disk accretion is in- 
stantaneously truncated at a time tj (since the disk size and 
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Fig. 19. — Evolution of star formation rate (relative to that at f = 0) and (B-V) optical colors of galaxies of different initial baryonic masses and gas 
fractions (as labeled), in high-resolution hydrodynamic simulations. Left: Evolution of an isolated quenched or "truncated" disk which is completely cut off from 
accretion/external gas supplies at ; = 0. Center: Evolution of a merger remnant after the final galaxy coalescence/starburst phase axt = 0. No feedback is included 
(i.e. the decrease in star formation rate derives entirely from gas exhaustion and shock-heating). Right: Evolution of a merger remnant, with feedback in the form 
of starburst-driven and quasar-driven winds. The mergers rapidly redden to the red sequence ((B — V) > 0.8) in < 1 Gyr, but an isolated disk (even with secular 
instabilities operating) remains blue. 

rable to that used here), and this is probably only relevant for 
the highest-redshift mergers. 

Two conclusions emerge from Figure [19] In the case of a 
truncated disk, the decline in star formation rate is gradual, 
so the (B - V) colors redden very slowly. Even in a truncated 
IO'^Mq disk after 3 Gyr, the galaxy colors are significantly 
bluer {{B-V) ^ 0.6) than those of a typical red sequence 
galaxy {{B-V) > 0.8). Furthermore, although our simula- 
tions allow for disk instabilities (see Figure 6 of Springel et 
al. 2005b), and do form spiral structure and even bars (seen 
in the small variations in star formation rate), this (even in the 
most massive, gas rich cases) does not consume sufficient gas 
to quench the disk. It is extremely difficult for secular mecha- 
nisms to exhaust all the gas, especially in the outer, low den- 
sity regions of disks, and only a small continued rate of star 
formation is necessary to keep the galaxies blue. 

Furtherm ore, it is o bserved locally and at redshifts up to 
z ^ I (Bell et all 120 04) that only a small fraction of galax- 
ies occupy the "green valley" between blue cloud and red se- 
quence. Assuming that ~ 1 /2 of ^ M, galaxies must cross the 
green valley since z = 1 (roughly what is expected from com- 
pariso n of the mass/luminosity functions, e.g. Marti n et"aLl 
|2007|), a slow reddening such as that of the truncated disk 
in Figure [19] would imply as many as ~ 1/4-1/3 of all 
^ galaxies sh ould occupy this region, compared to the 
< 10% observed (iBell et al.ll2004 . Simply put, this would 
eliminate or completely smooth out the observed strongly bi- 



baryonic mass should no longer evolve) - it will evolve as 
M oc (1 -I- [t-tf]/to)~^/^, where fo « r is a constant timescale 
which depends in detail on the gas fraction, gas mass profile, 
and time of star formation truncation {to = 0.72 r for a lO" Mq 
exponential disk with /gas = 0.4 truncated at z = 2). 

Figure [19] next shows the star formation rate and colors of 
merger remnants, after the merger itself. We consider two 
cases - first, with no feedback (i.e. no stellar winds, and no 
black hole accretion or feedback) included, and second, with 
a standard, observationally calibrated and relatively mild pre- 
scription for starburst-driven winds (with a wind outflow rate 
roughly half the star formation rate; see iCox et al.ll2007l) and 
BH accretion and feedback (such that the BHs self-regulate 
at masses appropriate for the observed Mbh - c relation; see 
iDi Matteo et all 12005 *). Stai- formation rapidly falls by or- 
ders of magnitude after the merger, even in the "no feedback" 
case, as the majority of the gas supply has been rapidly con- 
sumed in a central starburst and much of the remain ing gas 
shock-heated into an X-ray halo (ICox etal.ll2006ah . With 
feedback, the suppression is even more complete, as stellar 
and quasar-driven winds clean up the last remaining traces 
of star-forming gas. In both cases, the remnants redden ex- 
tremely rapidly, requiring less than one Gyr to reach the red 
sequence. We do caution that, in mergers of extremely gas- 
rich disks, f eedback may be necessary to redden so rapidly - 
ISpringel eF al. (2005a) showed this for 100% gas disks; but in 
any case the level of feedback required is reasonable (compa- 
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modal color-magnitude distribution (e.g. iBaldrv et al.|[2004l) . 
at least in moderately massive galaxies (our simulations are 
obviously not meant to be applied to e.g. dwarf satellites). In 
contrast, even gas-rich merger remnants with no feedback red- 
den rapidly to the red sequence, with a timescale for redden- 
ing of < 1 Gyr that is completely consistent with the observed 
color bimodality and small fraction of galaxies in the "green 
valley." 

4.2. "Maintenance ": How Is Later Star Formation 
Suppressed? 

It is apparent that merger remnants redden rapidly onto the 
red sequence. However, whether or not they can stay on the 
red sequence for significant periods of time is less certain. In 
other words, although mergers easily terminate star forma- 
tion, they will not remain as long-lived red galaxies unless 
they also maintain low levels of accretion and star formation. 

4.2.1. Is there a need to do so? The "No Feedback" Solution 

One possibility is that this maintenance is trivial. Roughly 
half the present mass density in red galaxies is built up since 
z = 1, and the typical host halo of a ~ M, red galaxy at z = 0.5 
will grow only by ~ 0.2 dex to z = 0. While this is still 
enough fractional growth (s^ 50%) to make the galaxy blue, 
if all the newly accreted gas were to cool immediately and 
form new stars, it is unlikely that this small amount of gas 
at the virial radius could cool and infall within the ^ 5 Gyr 
timescale to z = 0. Indeed, the "cooling flow problem" appears 
to be a prob lem for only the most massive cluster s at low red- 
shifts (e.g. Best et al. 2006; 'VikhH nin et al.ll2007l but see also 
[Chen et al. (2007)), which suggests that cooling flows are, in 
general, a late-forming phenomenon just now becoming rele- 
vant, and perhaps were never suppressed in the past. Further- 
more, many "central" galaxies do not actually reside at the ex- 
act center of their group or cluster potential (Mu lchaev et al.l 
12006; Jeltema et al. 2006), as is naively assumed in most an- 
alytic models, which makes the formation of coolin g flows 
less efficient. Recent hi gh-resolution simulations ( Naab et al.l 
120071: iKeres et al.ll2007h do suggest that, without any AGN or 
stellar feedback, the combination of virial shocks, compres- 
sion, and kinematic heating by clumpy accretion flows can 
prevent substantial cooling at z < 1, and that simple gas ex- 
haustion via star formation can quickly eliminate most of the 
low cooling-time gas. 

However, this is not entirely satisfactory, at least in the 
simplest sense. First, cooling flows still do appear to be a 
problem in these systems - and the most m assive galaxies 
are almost uniforrnly red (iBaldrv et al.l l2006i) . they do not 
appear to b e recently accreting/ star-forming or "becoming" 
blue (but see lRaffertv et al.l2006 l, who reach the opposite con- 
clusion for BCGs with large cooling flows). Second, even 
in the moderate-mass halo case considered above, the free- 
fall time of the gas accreted since z = 0.5, ^ 2 Gyr, is suffi- 
ciently short that a cooling flow problem remains a possibil- 
ity. The problem also becomes more severe at high redshifts, 
where cooling rates can be a factor ^100 higher than at z = 
(scaling oc n oc (1 +z)^). Finally, essentially all implementa- 
tions of galaxy formation models which attempt to account 
for gas accretion and cooling with various prescriptions have 
found that feedback of some kind is necessary to prevent new 
accretion in massive galaxie s (e.g. Birnboim & Dekel 2003; 
Binnev 2004'; 'Granato et al. 2004; Scannapieco & Oh 2004; 
Keres etal, 2005; Monaco & Fontanot 2005' CrotonetaLi 
2006l:lDekel & Birnboimll2006l:ICattaneo et al.ll2006l) . 



4.2.2. Can Quasar/Starburst Feedback Completely Suppress Future 

Cooling ? 

It is therefore natural to examine the feedback effects in- 
volved in (or stemming from) major mergers. We identify 
four primary feedback mechanisms: 

(1) "Kinematic" Feedback: Mergers themselves stir large 
quantities of gas, allowing relatively hot and cold gas from 
the inner and outer regions of a shared halo to mix, and 
generally increasing the cooling time significantly and dis- 
rupting any cooling flow s ongoing or in formation. This is 
seen bofli in simu lations jNaab et al.ll2007l: iKeres et al.ll200'7t 
[Cox et all 12006a') and X -ray observation s of galaxy groups 
(Jeltema et al. 2006; Vikhlinin et al.l2007l) . Furthermore, tidal 
shocks in the merger itself heat a significant quantity of gas to 
temperatures well above those that can efficiently cool in a 
Hubble time (the reason for the low star formation rates at 
late times in Figure[T9] even in the "no feedback" case). 

(2) Starburst- Driven Winds : It is known from local mea- 
surements (e.g. Kennicutt 1998) and also suggested in high- 
redshift studies (Erb et al. 2006) that a high surface den- 
sity of star formation inevitably results in strong galactic 
winds. Presumably driven by a combination of young stel- 
lar winds and supernovae, the energy coupling efficiency 
appears to be high (order unity), and simulations demon- 
strate that the observations are well-reproduced for reason- 
able, theoretically expected mass-loading efficiencies (r/ ^ 
0.5, where Mwind = riM^,, with possible mass depen dence from 
a momentum-based escape ap proximation; see ICox et akl 
120071; lOppenheimer & Pavel l2006l) . These wifl act through- 
out a merger, and are a powerful integrated source of feed- 
bac k, although not a s impulsive as quasar-driven outflows 
(e.g. lLidz et al]l2007h . 

(3) "Quasar" Feedback: Quasars ar e known to often ex- 
hibit strong outflows (for a review, see Veilleux et al]l2005l) 
and to have a large effect on the ionization and tempera- 
ture state of the inner regions of their host galaxies (e.g. 
iLaor et al.lll997l;lKrongold erani2007l;lRupke et al.ll^^ and 
references therein). At the brightest luminosities, a large frac- 
tion (~ 40%) of sightlines to quasars see highly energetic 
(> several lO^kms"') broad absorption line (BAL) outflows, 
and it is likely that all bright quasars exhibit some such BALs 
(alth ough they may not be visible owing to geometric ef- 
fects; iReichard et al.ll200l lEiviil2000l; iGallagher et alll2006l; 
iPriddev et al.l 120071) . These feedback mechanisms must be 
able to have a dynamical effect in some sense on the host, in 
order to suppress accretion onto the central BH once it reaches 
the limit of the Mbh - cr relation. Recall, if only ~ 0. 1 % of the 
initial galaxy gas mass were to survive a merger and make its 
way to the center of the galaxy, it would (without BH feed- 
back) drive the central BH off the Mbh - c relation by more 
than the observed scatter (while having almost no effect on 
a). 

What is less certain is the effect such feedback has on the 
largest galactic scales. In most models, it is inevitable that 
the small-scale wind, heating, or pressurization required to 
halt accretion and produce the Mbh - c relation will indeed 
generate a galactic outflow, and recent high-resolution, self- 
consistent simulations of quasar feedback and accretion disk 
winds imply the formation of powerful kinematic outflows 
that will couple on larger scales (Proga 2007). Indeed, an 
increasing number of bright quasar hosts have now been ob- 
served in which jets or wind s appear to be strongly impacting 
the host galaxy and halo gas (IZirm et al.l2005l : lNesvadba et al.l 
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120061 120071: iReuland et ani2007h . or in which BAL quasar 
winds are ent raining gas at ^ 10.^ .~ 10'*kms~' velocities o n 
~kpc scales (Ide Kool et alJ l200lll200l iGabel et al.ll2006h . 
Indeed, high velocity winds at (or beyond) galactic scales ap- 
pear to be ubiquitous in post-starburst galaxies at moderate 
redshifts, and trace a continuum in outflow properties with 
bright quasars (typically with low-luminosity AGN consis- 
tent with fadin g from a recent pe ak of starburst a nd quasar 
activity; Tremo nti etaP 12007: Gan guly et ani2007h . The in- 
tegrated energy in this feedback is, from simple energ etic ar- 
gume nts, comparable to that from stellar winds ( Lidz et alj 
l2007h . However, the timescale is much shorter - the BH 
gains most of its mass (releases most of its energy) in less 
than a Salpeter time ~ 10^'^ yr, whereas most of the stars are 
typically formed over a timescale > 10'^ yr Therefore, even 
in the most conservative models, the power in such quasar- 
driven winds is ^ 10-100 times greater than that in typical 
starburst-driven winds. 

As a result of the short timescales associated with this 
process, and because the energy or momentum is injected 
on scales small compared with those of entire galaxies, the 
impact of quasar feedback is e xplosive in nature. In deed, 
iHopkins & Hernquis j (l2006h and lHopkins et all (l2006b() have 
demonstrated that the outflows in the simulations caused by 
this phenomenon are well-approximated by a generalized 
blast-wave solution. 

(4) "Radio-Mode" Feedback: As coined by lCroton et al.l 
(|2006|) . this refers to a maintenance mode of feedback, includ- 
ing e.g. the inflation of r adio bubbles in clusters at relatively 
low accretion rates (e.g. Fabian et al. 2006: Du nn & FabianI 
|2006: Allen et al. 2006b: Sanders & Fabian 2007^ and refer- 
ences therein), but also the driving of weak winds from ra- 
diatively inefficient accretion flows (e.g. lNarayan & Yilll994l) 
and X-ray heating of nearby gas. Essentially, this is a blanket 
term for all feedback mechanisms which depend on a mas- 
sive BH at relatively low levels of activity (low Eddington ra- 
tios), in which most massive BHs spend most of their lifetimes 
(typically since z ^ 2). It explicitly does not include "quasar- 
mode" feedback, representative of the high-Eddington ratio, 
high-power output effects described above (and it also does 
not exclude radio jets as a feedback mechanism in the "quasar- 
mode"). 

This mode of feedback operates over long timescales fn) 
and would not occur under the high-Eddington ratio condi- 
tions of mergers. Nevertheless, we include it here because it 
is linked to mergers in a critical way: almost universally, the 
forms of "maintenance" fe edback require the presen ce of a 
relatively massive BH (e .g.lDekel & Birnboim 2006: Binneyl 
I2004t ICroton et al.ll2006t [Siiacki et al.' 2007). A massive BH 
empirically requires a massive spheroid, which requires a ma- 
jor merger. In other words, there is no "radio-mode "feedback 
without major mergers. Indeed, a comparison of the local 
mass density of supermassive black holes with the luminos- 
ity density of quasar s integrated over redshift dSoltanlll982l: 
iHopkins et al.ll20()79) indicates that the mass growth thro ugh 
the radio mode must be negligible (Hopkins et al.ll20b6dl) . in 
agreement with prelimina ry results from cos mologica l simu- 
lation s of these processes dSiiacki et al.l2007tiDi Matteo et alj 
l2007h . 

Although the distinctions between these modes of feed- 
back, and ultimately the detailed identification of the drivers 
of each are of great importance, these are questions outside 
the scope of this paper A detailed comparison, for example, 
of the effects of different modes of feedback on the IGM and 
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Fig. 20. — Top: Cooling time as a function of radius for typical relaxed 
SPH simulation merger remnants embedded in massive gaseous halos. We 
show both the absolute value of the cooling time (left) and the cooling time 
relative to the local free-fall time (right). Solid (dot-dashed) lines show 
the mass-weighted mean value at each r for galaxies with stellar masses 
^ 1O'-M0 (~ 3 X 10"'Mq). Different colors correspond to different ini- 
tial halo gas profiles: black assumes a relatively low total halo gas mass, red 
and blue a factor of several higher gas mass (with a pre-merger isothermal or 
pressure-supported temperature profile, respectively). For one case we show 
contours of the full gas distribution at 50, 5, 1 and 0.1% - most of the gas 
is close to the mean value. Vertical dotted lines show the virial radii of the 
halos of both masses at z = (black) and z = 2 (green). Lower: The integrated 
gas mass below a given cooling time (left) or cooling time relative to free-fall 
time (right). Dotted horizontal lines show the two galaxy stellar masses rep- 
resented here. Dashed line in left panels shows the Hubble time. Feedback 
from a major merger heats a large quantity of gas and establishes a hot or 
"quasi-static" halo. 

their observable signatures will be the topic of a future paper 
(Hopkins et al., in preparation). For now, we simply wish to 
examine whether a reasonable integrated effect of such feed- 
back could be to completely suppress future star formation in 
merger remnants. To do so, we return to the numerical simu- 
lations described above. 

First, we consider several merger simulations similar to 
those in Figure [191 with reasonable, observationally con- 
strained feedback prescriptions from both star formation and 
quasars. We embed the progenitor systems in large gaseous 
halos meant to represent the gas both in the host halo and sur- 
rounding it (which will be accreted in the future). Specifically, 
we consider a range of gas mass for the halo, from the total 
gas mass expected within 7?vii to several times this quantity 
(representative of that which will be accreted in the next few 
Gyr). We consider both an NEW profile for the gas and a uni- 
form density distribution, and calculate the initial gas temper- 
ature either assuming a uniform heating to the virial temper- 
ature or initial hydrostatic equilibrium. The halos are initial- 
ized appropriate for those at low redshifts, and the galaxie s are 
otherwise constructed identically to those described in § 14.11 
We evolve our merger simulations until they are relaxed (typ- 
ically '^2-3 Gyr, as before), and calculate the cooling time 
(including metal-line cooling, following Cox et al. 2006a) for 
all remaining gas particles at the end of the simulation. 

Figure|20]shows the results for several representative simu- 
lations. For clarity, we do not show the results of every sim- 
ulation, but note that the qualitative behavior is, in all cases 
of a given mass, quite similar, with properties such as the ini- 
tial gas density profile affecting only the details of the final 
gas profiles (not thek general behavior as a function of mass 
and/or radius). We plot the gas cooling time as a function of 
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radius - by this time, the gas has relaxed and there is a rea- 
sonably well-defined cooling radius inside which the gas will 
cool in <C tii. The actual mass contained therein is not neg- 
ligible - only ~ 5 - 10% can be added to the galaxy mass in 
a Hubble time. This is, in principle, sufficient to make the 
galaxy blue once again, however most of the cooling would 
happen at late times - where another small burst of feedback 
could re-heat the gas and prevent this scenario. For galaxies 
moving to the red sequence at z < 1 (~ 1/2 of present red 
galaxies), the suppression is even stronger - only < 1 % of the 
post-merger galaxy stellar mass can cool by z = 0, which is 
sufficiently small to ensure that the galaxy remains "red and 
dead." 

These simulations, however, neglect the dynamical nature 
of accretion onto the dark matter halo with cosmic time, and 
do not include the very large relative gas mass accreted onto 
the most massive, early-forming systems. For example, a > 
few X lO'^M© BH forming at z = 2 will live in a > 10'^ /;"' Mq 
halo, which will typically grow by more than an order of mag- 
nitude in mass to z = 0. In our complete SPH simulations, 
surveying this parameter space requires large boxes, external 
reservoirs of gas, inclusion of cosmological effects, and long 
runtimes tu), and is ultimately outside the scope of this 
paper. However, we can make some rough estimates of the 
qualitative effects from simple scaling arguments. 

Consider the feedback energy which couples to the galac- 
tic ISM during a merger (/Smerger)- The integrated feedback 
energy injected by the "quasar mode" over the course of the 
merger will be (given that most of the BH mass is gained in 
this phase) approximately E = rjerMBHC^, where ~ 0.1 is 
the radiative efficiency and rj is the feedback coupling effi- 
ciency (rj 0.05 in our simulations in order to yield the appro- 
priate normalization of the Mbh - relation). As noted above, 
the total feedback energy from star formation is of a compara- 
ble order (although it operates over a much larger timescale, 
so only some fraction will couple during the merger itself), 
so we can subsume it into this scaling (since ?/ is uncertain 
anyways, it can effectively include the maximal factor ^ 2 
addition from stellar winds and supernovae). The feedback 
from BH growth and at least the final, peak starburst phase 
will couple in a short time, ~ 10^ years (the timescale for 
the final e-foldings of BH growth), much shorter than the dy- 
namical time in the outer regions of the halo. 

Assuming, therefore, that the "merger feedback" (by which 
we mean the combined feedback from quasar, starburst, and 
kinematic effects - although the latter are energetically sub- 
dominant) creates a strong shock (true in nearly all of our 
simulations), the post-shock temperature inside the virial ra- 
dius of the halo will be approximately Tshock ~ ctcTyi^, where 
c is the halo concentration and a is a coefficient of order 
unity which depends in detail on th e halo gas profile, baryon 
fraction, and metallicity (we follow lOekel & Birn boim 2006, 
who adopt standard values for these quantities, and obtain 
a w 0.5). The resulting cooling time of the shocked gas near 
the virial radius is then (for the same parameters) roughly 
fcooi «^ 8.3A2;)(l+z)-^(7;hock/10^K)2fH, where A200 ~ 1 is 
the virial overdensity at the given redshift relative to a value 
of 200 (and we approximate the cool ing function around 
the te mperatures of interest following I Sutherland & Dopltal 
Il993h . (Of course, most of the gas relevant for cooling will be 
at smaller radii and A200 ^ 1 , but we simply wish to illustrate 
the relevant scalings.) 

If we were to rely on one feedback event alone to suppress 
all cooling until z = 0, we would require two basic criteria. 



First, this clearly requires that Tshock be sufficiently high such 
that fcooi > tniz), i.e. 7;hock > TLritH « 4 x 10^K[A2oo/h (1 + 
2)3] 1/2^ where fn is the fractional lookback time to redshift 
z. Second, the coupled feedback energy must be sufficient to 
heat all of the total z = halo gas content to these temper- 
atures - i.e. the total mass which can be shocked, Mshock = 
^ifip E^ei-gf^r /{3 /IkTcrn.n) (where ^ = 0.59 for pristine gas) 
must be equal to or greater than the total gas mass which will 
be accreted by z = and therefore which must be prevented 
from cooling. The first criterion is satisfied for all moder- 
ate halo masses of interest (although it may be that low-mass 
halos at high redshifts z > 2 have difficulty shocking to suf- 
ficiently high temperatures), and this is borne out by direct 
comparison with the post-shock temperatures in our simula- 
tions at all redshifts. It is at least likely that some of the sur- 
rounding gas will be shocked to very high temperatures - the 
more interesting question is how this mass compares to the to- 
tal mass that will be accreted and (potentially) otherwise cool 
by z = 0. 

Given our expectation for the average galaxy, and corre- 
sponding BH mass, in a halo of a given mass at some red- 
shift, Figure|2T|compares the mass that can be shocked (given 
the energetic criterion above) to that accreted by z = 0. In 
all cases, the feedback is able to shock-heat up to several 
times the initial galaxy mass, and we crudely expect the shock 
to propagate to several times the initial virial radius of the 
galaxy. However, the implications of this can be quite dif- 
ferent for halos of different masses. Low mass halos, even 
at z = 2, grow by a relatively small amount. For example, 
an average 10" /;~'Mq halo at z = 2 grows by a factor 5 
to z = (so the feedback from the merger need only shock 
several times the galaxy mass in external gas to prevent all 
future accretion), but an average 10^^ Mq halo at z = 2 
grows by a factor ^ 25 to z = 0. In small halos, then, feed- 
back from a merger, at least at redshifts z < 2, may be able 
to completely prevent future accretion, without the need to 
invoke any maintenance mode of feedback. In large halos, 
however, there is too much continued accretion and growth at 
low redshifts, and there is little chance that a single, merger- 
triggered burst of feedback can (alone) suppress all future 
growth. The division between the regimes appears to be at 
Mhd\o ^ 10'^ - 1O'-'/!"'M0, interesti ngly s imilar to the tra- 
ditional halo quenching mass (see § 14.2.31 below). We note 
that this analysis can be repeated in terms o f the post-shock 
entropy, following'S cannapieco & Ohl (|2004|) . which yields a 
nearly identical result. 

At the highest redshifts z > 2, it is also difficult for a single 
event to suppress the cooling of all gas which will be accreted 
by these halos, especially for systems which are already mas- 
sive at these redshifts (and therefore likely to form the most 
massive clusters at z = 0). A ~ 10'^/i~'Mo halo at z = 4, for 
example, is likely to grow to a ~ 10'^ Mq cluster by z = 0, 
so the baryon content contributing to the merger-driven feed- 
back event at these redshifts is negligible compared to that 
which will be accreted at later times. We also caution that the 
feedback from a merger may not be as efficient as we have 
assumed in this analysis. Although we adopted a relatively 
conservative total "stellar-nquasar" feedback energy input, it 
is not entirely clear how successful such feedback is at cou- 
pling to gas on large scales. Perhaps more important, the sim- 
ple scalings above ignore the possibility that cooling instabil- 
ities might occur within the post-shock compression, or that 
cold clumps might be able to self-shield against a propagat- 
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Fig. 21. — Left: Total gas mass (relative to the total gas mass of the halo 
at z = 0) which can be shocked by a merger-induced feedback-driven out- 
flow/blastwave at the given redshift and host halo mass above temperatures 
for which the cooling time becomes longer than the Hubble time (i.e. total 
fraction of the z = halo baryon content for which cooling can be completely 
suppressed until z = by a typical merger at the given redshift). Right: To- 
tal gas mass (relative to the total gas mass of the halo at the given redshift) 
which can b e shocked to the critical shock stability temperature (following 
IKeres et al] {2005); Dekel & Birnboim 1 2006)) above which the cooling time 
is longer than the (instantaneous local) gas compression or free-fall time and 
a quasi-static halo is established. Feedback from a major merger alone can 
quench all future accretion in halos below the traditional "hot halo" mass 
threshold at moderate redshifts (z < 2), and can easily estabhsh a "hot halo" 
within the virial radius at all redshifts. 

ing shock, leaving most of the mass which would be accreted 
unaffected. 

In particular, if gas accretion occurs preferentially along 
filamentary structures, it may be difficult for feedback to di- 
rectly couple to most of the gas in the filament. A more de- 
tailed calculation of these effects will, unfortunately, require 
better knowledge of the actual drivers of feedback, as well 
as high-resolution simulations which can self-consistently re- 
solve phase structure and shocks in the IGM gas. For now, we 
would more cautiously describe our calculations as estimates 
of what feedback from a major merger could do to suppress 
cooling. Even in this case, however, both our SPH simulations 
and simple scaling arguments suggest that the most massive 
systems, especially if their mergers occur early at z > 2, can- 
not be quenched just by the energy injection from a single 
feedback event. 

4.2.3. A "Mixed" Solution: Hot Halos from Quasar/Starburst 
Feedback 

Given the uncertainties and limits, in the most massive sys- 
tems, on the efficiency of short-term feedback from a major 
merger, we propose a mixed solution. Halos more massive 
than (roughly) ^ 1O'^/!~'M0 have characteristic gas cooling 
timescales longer than the dynamical or free-fall timescale, 
and are so described as being in the so-called "quasi-static" 
or hot halo regime. In most models, "radio mode" feedback, 
i.e. some form of feedback from low accretion rate activity in a 
massive central BH provides the small additional heating term 
needed to maintain a pressure-supported hydrostatic equilib- 
rium structure at all radii, preventing new gas from cooling 
onto the central galaxy. The additional heating term, for our 
purposes, does not even necessarily need to come from a cen- 
tral BH - it could owe to kinematic heating or other effects, 
so long as i t maintains the hot halo - although energetic argu- 
ments (e. g. iBenson et al.l[2003;) a nd high-resolution observa- 
tions (e.g' lBatcheldor et ani2007l) favor an AGN origin. 

A significant pro blem with these models, however, as we 
have seen in the Cro ton et al.l ([2006) example in § 13.21 is that 
they are unable to produce sufficient numbers of red central 
galaxies in relatively low mass halos, and the red fraction 
does not depend on stellar mass as is observed. In other 



words, some process, with a dependence on galaxy mass, is 
required to assist the quenching of lower-halo mass systems. 
One might attempt to address this by adding a strong secular 
quenching mechanism, but we have shown by including this 
in the Bower et al. (2006) models that this fares little better at 
matching the bivariate red fraction as a function of halo and 
galaxy stellar mass, and that it conflicts with constraints on 
pseudobulge populations. 

However, we have just shown that feedback from a major 
merger can shock-heat sufficient surrounding gas to quench 
systems below the traditional hot halo mass threshold (~ 
10'^- IO'-'Mo) for substantial periods of time. We therefore 
propose that traditional modes of quenching and feedback in 
hot halos remain the key to suppressing star formation in mas- 
sive systems, but that these are supplemented by mergers, 
which can effectively quench star formation in lower-mass 
systems before these cross the hot halo threshold. In fact, the 
major merger needs to suppress star formation in low mass 
systems only until they would naturally develop hot halos - 
often much less than a Hubble time. For example, a typical 
^ 10'' /z~'M0 halo merging at z ~ 4 need only be quenched 
by merger feedback until z ~ 2 (« l.SGyr), when it will be 
sufficiently massive to enter the traditional hot halo regime. 
Once a hot halo is developed, the merger remnant already, by 
definition, has the means to maintain that halo and supplement 
it with feedback - namely, a relatively massive spheroid and 
BH which will be accreting at low rates (i.e. the ideal seed for 
"radio-mode" feedback). 

More conservatively, the "merger feedback" does not even 
need to completely suppress cooling/accretion in these low 
mass systems. If the hot halo is an effective means of quench- 
ing, then mergers only have to create hot halos. In fact, the 
traditional hot halo is generated by an accretion shock in mas- 
sive systems, and does not occur in low-m ass systems because 
the co nditions do not set up such a shock dPekel & Birnboiml 
l2006h . It is a small extension, then, to suppose that the 
strong shocks from merger-induced (quasar and starburst- 
driven) feedback, which are powerful even in low-mass sys- 
tems, might accomplish this even when accretion shocks do 
not. Indeed, in Figure |20] we show the cooling time rela- 
tive to the free fall time (fff) for the gas in our SPH merger 
remnants, and the amount of gas mass raised by the feedback 
coupling above a given fcooi/fff- Regardless of the mass of the 
systems or absolute values of the cooling times, the gas out 
to many times the virial radius is almost uniformly raised to 
fcooi/fff ^ 1, the traditional criterion for a hot halo. 

In other words, relatively low-mass halos (which would 
otherwise rapidly cool) require some event to enable their 
quenching and transition to a stable hot accretion mode (i.e. 
suppression of future cooling). We demonstrated in § 14.2.21 
that feedback associated with a major merger can easily ac- 
complish this (although we note that feedback may be inef- 
ficient in extremely low mass halos, < lO'^M©, which are 
not important for our conclusions). Once a halo grows to 
large masses, however, any single quasar or stellar feedback 
event (or any baryonic feedback event, given the relative mass 
growth involved) is probably insufficient to singlehandedly 
heat the (very large) quantities of gas involved to a temper- 
ature so that the cooling time is longer than a Hubble time. 
However, in this regime, massive halos are already heating 
most of the gas via accretion shocks. Some mechanism (such 
as a major merger) is still needed to exhaust the gas in the cen- 
tral galaxy, and additional mechanisms (such as radio-mode 
feedback; requiring a massive spheroid and black hole in the 
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remnant) may be needed to account for a small additional en- 
ergy input or mixing term in the center of the halo (in order to 
prevent the formation of cooling flows at late times), but the 
bulk of the energetic input needed to maintain a hot system is 
already in place. 

We specifically check this scenario by revisiting Figure 1211 
and considering, instead of the amount of gas which can 
be heated to temperatures above which the cooling time is 
much longer than a Hubble time, the amount of gas which 
can be shocked to temperatures which are above the critical 
shock stability threshold; i.e. for which the cooling time is 
lon ger than the free-fall or dyn amical time of the gas. Follow- 
ing lDekel & BirnboimI (l2006l) . we estimate this critical tem- 
perature for gas near the virial radius of the halos of inter- 
est, and obtain the (halo-mass independent) threshold Tan ~ 

4 X 10^KA2qq(1 +2)^/"*. Interestingly, the ratio of BH to host 
mass (and therefore the relative amount of feedback energy 
coupled to the gas) appears to scale with redshift in a roughly 
similar manner (see Paper I and Hopkins et al. 2007c), yield- 
ing a nearly redshift-independent ratio of mass which can 
be shocked by merger feedback to that inside the virial ra- 
dius at each epoch. In other words, at all redshifts, feedback 
from quasar and/or starburst activity associated with a ma- 
jor merger is sufficient to shock the entire gas content within 
the virial radius (or even to several times the virial radius) to 
this critical temperature, for halo masses M^mo ^ IO'-'Mq. At 
larger halo masses, systems will already have naturally devel- 
oped hot halos owing to accretion shocks, so it does not matter 
whether or not the feedback energy can shock the systems into 
the hot halo mode (although the merger-driven exhaustion and 
feedback may still be critical to ceasing star formation and 
making the system red). And once the hot halo mode is es- 
tablished in the inner radii inside /?vir, it does not ultimately 
matter how far the hot halo extends beyond the virial radius 
(or how much of the mass to be later accreted is affected) - 
the hot halo sets up a quasi-static, pressure supported equilib- 
rium against which newly accreted gas will shock and add to 
at large radii (regardless of its mass). 

A more detailed study of these hot halos from ma- 
jor merger-driven quasar and starburst feedback in cosmo- 
logical simulations is an important topic of future work. 
However, it is ultimately a relatively small variation on 
the traditional principle which has been recognized for 
many years (see Rees & Ostrik e^ll977l: lN orman &Silkl[T97l 
iBlume nthal et al. 1984). We have further shown that it is not 
only possible, but quite easily accomplished from moderate 
feedback prescriptions. Quenching can therefore be accom- 
plished in the "traditional" context of hot halos supplemented 
by feedback from a massive BH, but allowing for feedback 
from black hole growth and star formation in a major merger, 
in a halo of any mass, to create a hot halo environment. 

5. DISCUSSION 

We have developed and tested a simple but physically- 
motivated model in order to study the cosmological role of 
mergers in the formation and quenching of red, early-type 
galaxies. By combining theoretically well-constrained halo 
and subhalo mass functions as a function of redshift and envi- 
ronment with empirical halo occupation models, we can pre- 
dict the distribution of mergers as a function of redshift, envi- 
ronment, and physical galaxy properties. In Paper I, we dis- 
cuss this methodology in detail, and show that it accurately 
reproduces a variety of observations over a wide range in 
redshifts, including observed merger mass functions; merger 



fractions as a function of galaxy mass, halo mass, and red- 
shift; the mass flux/mass density in mergers; the large-scale 
clustering/bias of merger populations; and the small-scale en- 
vironments of mergers. The primary advantage of this model 
is that it allows us to study and make a priori predictions for 
the effects of mergers without many of the uncertainties or 
degeneracies inherent in present cosmological simulations or 
semi-analytic models. 

For example, cosmological simulations still lack the resolu- 
tion to model the processes of internal galactic kinematics in 
mergers, black hole accretion/feedback, and disk formation. 
Although progress is being made studying these processes via 
"zoom-in" simulations, it is not meaningful to speak of gas- 
rich, spheroid forming mergers in cosmological populations 
if a cosmological box does not contain the appropriate, rep- 
resentative population of accurately formed disk galaxies (the 
progenitors in these mergers) in the first place. 

Although semi-analytic models avoid some of these diffi- 
culties, they require making a number of assumptions regard- 
ing models or physics that we are not attempting to test in this 
paper, including e.g. disk formation, star formation efficiency 
in disks, disk instabilities, minor mergers, satellite disruption, 
reddening of satellite galaxies, and the exact physical mecha- 
nisms of feedback. These assumptions introduce uncertainties 
in the model and, more importantly, obscure the key physical 
elements being tested. 

Our adopted model, in contrast, bypasses these (unneces- 
sary for our purposes) assumptions and uncertainties, and in- 
stead empirically adopts the relevant consequences of all these 
physical processes - namely what kinds of galaxies are merg- 
ing at a given place and time. We can then more directly ask 
the question we wish to answer: how do mergers contribute 
to the formation and/or quenching of massive red galaxies? 

We find that the simple assumption that star formation is 
quenched after a gas-rich, spheroid-forming major merger 
(by any mechanism) naturally predicts the turnover in the 
galaxy mass-halo mass relation at ~ - i.e. the fundamen- 
tal turnover in the efficiency of star formation and incorpora- 
tion of baryons in galaxies, at the observed scale and with- 
out any parameters tuned to this value. The physical scale 
^ reflects the point where major, galaxy-galaxy merg- 
ers first become efficient. At lower masses, major mergers 
are rare - this is true bot h of halo-halo major mergers (e.g. 
Ivan den Bosch et al.l2005h and galaxy-galaxy mergers (which 
are further suppressed at these masses because of the relative 
scalings of orbital velocities and internal galaxy velocities - 
i.e. two such galaxies are likely to interact as field flyby or 
satellite-satellite systems with relatively high orbital veloci- 
ties that do not efficiently merge). 

Systems therefore generally grow uninterrupted, potentially 
building relatively low-mass pseudobulges (< 10'"Mq) via 
disk/bar instabilities or minor mergers, until they get to ^ . 
By these masses, the probability of the halo merging with 
a major companion reaches of order unity, and the veloc- 
ity scalings are such that the two galaxies (once the halos 
have merged) will merge efficiently (fmerger <C fe)- The sys- 
tems can then grow via subsequent (dry) mergers, but this 
is a relatively inefficient channel (i.e. mass growth is slow). 
Because their star formation is quenched (and therefore no 
longer keeping pace with their host halo growth), mergers 
themselves also rapidly become less efficient (i.e. the system 
mass becomes low relative to the host halo mass, increasing 
the merger timescales). 

In addition, our model naturally predicts the observed mass 
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functions and mass density of red galaxies as a function of 
redshift, the formation times of early-type galaxies as a func- 
tion of mass, the fraction of quenched galaxies as a func- 
tion of galaxy and halo mass, environment, and redshift, and 
the distribution/dichotomy of kinematics in massive ellipti- 
cals. Each of these predictions agrees well with observa- 
tions over the entire observed range of galaxy masses and 
redshifts. As demonstrated in Paper I, our model also agrees 
well with observed merger rates and fractions as a function 
of galaxy mass and halo mass at all observed redshifts. To- 
gether with the agreement between our model and the ob- 
served mass functions and mass density of red galaxies, this 
illustrates that there are, in fact, sufficient numbers of merg- 
ers (both in theory and observed) to produce the entire mas- 
sive sphero id population at all observed redshifts (see also 
iHopkins e t al. 2007a). Also, unlike commonly adopted mod- 
els in which quenching is regulated purely by halo mass, we 
have not adjusted or tuned any parameter to give the desired 
results. Indeed, there is not even an obvious parameter which 
can be tuned to give the turnover in the galaxy mass-halo 
mass relation at the appropriate location (since it appears not 
to depend on our calculation of the merger timescale). To 
the extent that mergers can supplement quenching, then, this 
suggests that it is not necessarily problematic that theoret ical 
calculations (Bi rnboim & Dekelll200l iKeres et al.ll2005h do 
not give exactly the same halo quenching threshold as semi- 
analytic models subsequently tuned to fit the observations, as 
has been noticed in several works (e.g. ICroton et al.l 120061; 
ICattaneo et ani2006l) . 

Although these predictions are suggestive, recent semi- 
analytic models have demonstrated that many of them are 
non-unique. A variety of different quenching implementa- 
tions and feedback effects in these models have been shown 
to successfully reproduce e.g. low-redshift mass functions, 
color-magnitude diagrams, and mean red fractions. We there- 
fore investigate the robust, observable differences between 
three broad classes of models for quenching. 

First, our adopted merger-induced quenching model, in 
which some mechanism enables merger remnants to remain 
quenched. Second, a halo quenching model, in which quench- 
ing is primarily determined by a simple (albeit potentially 
redshift-dependent) halo mass threshold (regardless of merger 
history or morphology) - i.e. one in which some mecha- 
nism enables any system to remain quenched if and only if 
it develops a "hot halo." Third, a secular model, in which 
color (and/or morphological) transformation is driven solely 
by galaxy structure (essentially baryonic galaxy mass), owing 
to e.g. disk or bar instabilities (or other non-merger related 
mechanisms). Regardless of the exact details of their quench- 
ing prescriptions (and other assumptions), most present semi- 
analytic models can clearly be identified with one of these 
three classes of models, based on which criterion effectively 
dominates quenching (e.g. galaxy merger history, halo mass, 
or disk mass), and we demonstrate that the key qualitative 
predictions of each class will remain true. Note that we are 
explicitly referring to the quenching of central halo galaxies 
(the great majority of > 1O"*M0 galaxies), as the reddening 
of satellites is almost certainly affected by other processes 
(such as their initial accretion, ram pressure stripping, or ha- 
rassment). 

We show that these models make a number of robust, 
unique predictions with respect to several observables, includ- 
ing: 

(1) Bivariate Red Fractions: Observational measurements 



of the red fractions of galaxies in groups can now break the 
degeneracy between the fraction of quenched systems as a 
function of galaxy mass (which all these models successfully 
reproduce) and halo mass. The observations show several 
important qualitative trends in the fraction of quenched, cen- 
tral halo galaxi es as a bivariate function of galaxy stellar and 
halo mass (e.g. IWeinmann et al.ll2006ah . These include: (1) 
a strong dependence of red fraction on halo mass, (2) some 
(weaker) residual dependence on galaxy mass/luminosity, (3) 
a lack of any sharp characteristic scale in Mhaio, (4) a relatively 
high red fraction (/red ^ 0.5) for the most massive/luminous 
systems even at relatively low halo masses (Mhaio ^ lO'^Mo), 
and (5) a similar, relatively high red fraction (/led ^ 0.5) 
for the least massive/luminous systems at high halo masses 
(A^haio ^ lO'^Mo). The fundamental difference between the 
classes of models we consider is directly reflected in this pre- 
dicted bivariate red fraction (where we refer specifically to 
central galaxies, as satellites may be affected by other pro- 
cesses as indicated above). 

In halo models, the red fraction is essentially a step func- 
tion in halo mass with a sharp transition from low red frac- 
tions to /red ~ 1 around the critical quenching mass, and little 
residual dependence on galaxy properties. In secular mod- 
els, the red fraction is just a function of galaxy mass, with 
little (or even inverse, if quenching becomes harder to main- 
tain in high mass halos) correlation with halo mass. Mergers, 
however, depend on both galaxy and halo mass, with larger 
galaxies at a given halo mass merging more efficiently (and 
being more likely to have already undergone a major merger), 
while larger halos are more evolved and more likely to have 
accreted a major companion as fuel for a major merger. Con- 
sequently, the red fraction is an increasing function of halo 
mass, but with an additional (weaker) dependence on galaxy 
mass, and grows smoothly (i.e. without a single, sharp char- 
acteristic scale) to higher masses. A significant dependence 
on halo mass is maintained, but there is still a large red frac- 
tion for the massive galaxies (even in relatively low-mass ha- 
los). More detailed observations are needed to quantify this 
in greater detail, but only the merger model appears to match 
the qualitative trends observed. 

(2) High-Redshift Passive Galaxies: A relatively large pop- 
ulation of massive, red galaxies exists at even high redshifts 
z > 3. Although at high redshifts most (simply identified) 
"red" galaxies are dusty star-forming systems, there is a sig- 
nificant population which are truly "red and dead," spectro- 
scopically confi rmed passively evo l ving, low star-form ation 
rate spheroids (lLabb6 et al.1 l2005t iKriek et al.l l2006l) . In 
contrast, semi-analytic halo quenching models are generally 
forced to assume that some process (e.g. accretion in filaments 
or cold clumps) at these redshifts raises the mass threshold 
for quenching, and as a result the predicted density of passive 
galaxies drops rapidly at z > 2. We note that this is not a state- 
ment that "hot halos" cannot or do not for m at these redshift s 
(simulations, in fact, suggest that they do; iKeres et ani2005h . 
nor that such models do not predict a sufficient density of all 
massive galaxies at these redshifts. However, in the naive im- 
plementation (in which the quenching is strongly dominated 
by a simple halo mass threshold), one cannot simultaneously 
form massive galaxies (and predict a sufficiently high global 
star formation rate density) at high redshift and quench them. 
In order to match both the observed density of star-forming 
and passive massive galaxies, some mechanism is required 
which can explain the quenching of some, but not all, systems 
in massive halos at high redshifts. 
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Mergers, on the other hand, proceed efficiently in mas- 
sive halos at high redshifts, predicting a significant density 
of quenched, passively evolving systems even at z > 3, in 
good agreement with the observations. A secular model can 
also explain the density of passive systems at these redshifts, 
since, by definition, the existence of such massive galaxies 
in the first place guarantees that a large fraction will be red 
(since the red fraction is a pure function of galaxy mass in 
this model). However, the secular model encounters a differ- 
ent conflict at high redshift. 

(3) Buildup of the Color-Density Relation: The color- 
density relation appears to weaken with redshift, flattening 
in intermediate density environments until z ~ 1.5, where it 
appears that there i s no measurable color-density relation in 
field environment s dNuiiten et al.l 120051 ICooper et al.ll2006al : 
iGerke et"ani2007h . Even at high redshifts z ^ 3, however, 
there is still a significant color-density relation (Oua dri et alj 
|2007|) - it is simply that the relation becomes significant only 
in more extreme (proto-cluster, for example) environments. In 
other words, a large population of quenched galaxies emerges 
rapidly at early times in the most massive environments, and 
then subsequently builds up in more moderate environments 
at lower redshifts. 

A halo quenching model predicts something similar to the 
low-redshift evolution in these trends (although with diffi- 
culty in producing quenched systems in all but the most truly 
extreme environments at high redshift, as described above). 
At higher redshift, systems above the halo threshold quench- 
ing mass represent progressively more extreme environments, 
and if this effective mass threshold increases with redshift, 
the trend is more pronounced. The red fraction is still nearly 
a step-function at each redshift, but with a shifting relative 
scale. 

A merger model also predicts a trend qualitatively similar 
to that observed. The most dense environments undergo their 
epoch of major mergers more rapidly than less dense envi- 
ronments (equivalently, more massive present environments 
passed through their small group stage at earlier times), al- 
though the red fraction in halos of all masses decreases with 
redshift (as there is less time for mergers to operate). By 
z ^ 1.5, typical field environments have uniformly low red 
fractions, and no significant measurable color-magnitude re- 
lation is expected. The location of the buildup of quenched 
galaxies shifts to denser environments, similar to the observed 
trend. 

A secular model, in contrast, predicts almost no evolution 
in the trend of red fraction with halo mass as a function of 
redshift (as a consequence of there being relatively little evo- 
lution in the average mass of a star-forming galaxy hosted by 
a given halo mass). If one allows for high-redshift disks being 
more compact, this increases their inferred instability, yield- 
ing opposite evolution in the red fraction versus halo mass 
to that observed (i.e. increasing quenched fractions with red- 
shift at fixed Mhaio)- As a consequence, although there may 
be some artificial evolution with redshift in the red fraction 
as a function of environment (as the same halo mass corre- 
sponds to different environments), there is no significant true 
evolution. Furthermore, in a secular model, the halo masses 
corresponding to field environments do not trend towards uni- 
formly low red fractions by z > 1 .5 - i.e. there is little "smear- 
ing out" of the color-density relations at high redshift. Future 
observations are needed to make these comparisons formal, 
but quantifying the evolution with redshift in the red fraction 
as a function of host halo mass (from large samples which can 



isolate groups and group central galaxies, and span a wide 
range of environments) will be a powerful discriminant be- 
tween these models. 

(4) Spheroid Kinematics (Dichotomy of Elliptical Galax- 
ies): Numerical simulations and observations of merger rem- 
nants and elliptical kinematics demonstrate that gas-rich ma- 
jor mergers (i.e. those involving disks, even with low gas frac- 
tions /gas ^ 0. 1) generally produce typical low-mass (< a few 
L*) ellipticals with central cusps, disky isophotes, and signif- 
icant rotation, while subsequent gas-poor spheroid-spheroid 
mergers produce typical high-mass ellipticals with central 
cores, boxy isophotes, and little rotation. There is a well- 
defined transition between the two classes of spheroids, at a 
mass ~ 2-3 X 10" Mq for each of these criteria. 

A merger model naturally predicts this transition point; at 
lower mass, most spheroids (i.e. merger remnants) have expe- 
rienced only their initial, disk-disk spheroid-forming merger 
or (in some cases) one additional, gas-rich, disk-spheroid ma- 
jor merger. At higher mass, most systems have undergone 
an additional, subsequent spheroid-spheroid major merger. If 
the major merger is associated with quenching, the low-mass 
disk mergers are guaranteed to be gas-rich, and the spheroid- 
spheroid mergers at high masses are guaranteed to be gas- 
poor, matching the observed trends and transition point in 
each of the cusp/core, disky/boxy, rapid/slow rotation crite- 
ria. 

In a halo quenching model, however, many systems un- 
dergo their first major merger somewhat before their host 
halos cross the quenching mass threshold, and therefore re- 
accrete significant disks. Their subsequent mergers are not 
gas-poor, spheroid-spheroid any longer, but gas-rich, disk 
mergers. As a result, the predicted transition mass between 
disky and boxy ellipticals is increased by an order of magni- 
tude (only the most massive cD galaxies cross the quenching 
mass threshold early enough to have had multiple subsequent 
major mergers since that time), in contradiction to the obser- 
vations. 

A secular model suffers from the opposite problem. In 
order for secular mechanisms to dominate quenching, they 
must act before major mergers transform the system to a 
spheroid - i.e. systems must (by definition in such a model) 
predominantly quench before they undergo their first merger. 
A large fraction of even these first mergers, then, are gas- 
poor, spheroid-spheroid (or pseudobulge-pseudobulge) merg- 
ers. The predicted transition mass between disky and boxy 
ellipticals is therefore decreased by an order of magnitude, 
again in contradiction to the observations. It appears that 
matching the observed transition in elliptical types (with- 
out violating the basic kinematic constraints from simula- 
tions and observations) fundamentally requires some quench- 
ing of massi ve spheroids/mer ger remnants (a conclusion also 
reached by Naa b et aT]|2006bl who begin from a halo quench- 
ing model). 

(5) Effects of Small- Scale Environment: It has also been 
suggested observationally that the red fraction (at fixed halo 
mass and galaxy mass) does not depend on large-scale en- 
vironment, but may depend on small-scale environment, in 
the sense that it incr eases with overdensities on small scales 
jBlanton et al.l 120061) . We caution that, at present, interpre- 
tation of these observations is difficult because they include 
both satellite and central galaxies. However, if the result is 
borne out for central galaxies alone (i.e. a central galaxy is 
more likely to be red, all else being equal, if it lives in a small- 
scale overdensity) via measurements of the cross-correlation 
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function for red, central galaxies and other galaxies, then this 
would also favor a merger model. We note that it is not neces- 
sary that merger remnants live in such overdensities long af- 
ter their mergers (as, by definition, the mergers will consume 
some of the very galaxies that define such an overdensity). 
However, if such a trend exists, it is difficult to explain in a 
pure halo quenching or secular model, as both mechanisms 
operate independent of neighboring galaxy populations. 

These consequences of merger-driven, halo mass-driven, 
and secular/disk instability-driven quenching models are ro- 
bust, and future observations should be able to break the de- 
generacies between the models. Although the quantitative de- 
tails may differ slightly in different implementations of the 
models, we have s hown that current state-of-the-art semi- 
analytic mo dels (e.g. lCroton et al.l l2006', Cattaneo et al."2006l: 
iBower et al. 2006, which include a number of other prescrip- 
tions and more detailed physical recipes than the toy mod- 
els described above) fundamentally yield predictions which 
are qualitatively identical to the behavior expected for the ba- 
sic classes of models described above. These behaviors are 
generic to any model in which these processes dominate the 
quenching of central galaxies, and the predictions shown are 
at least qualitatively robust regardless of "tuning" the models. 

That we have not tuned or adjusted our model to give a 
particular result should not, of course, be taken to mean that 
there are no uncertainties in our approach. However, we re- 
calculate all of our predictions adopting different estimates 
for the subhalo mass functions and halo occupation model 
(and its redshift evolution) and find this makes little differ- 
ence (a factor < 2) at all redshifts. The largest uncertainty 
comes from our calculation of merger timescales, where, at 
the highest redshifts (z > 3), merging via direct collisional 
processes may be more efficient than merging via dynami- 
cal friction, given the large physical densities. More detailed 
study in very high-resolution numerical simulations will be 
necessary to determine the effective breakdown between dif- 
ferent merger processes. Nevertheless, the difference in our 
predictions at these redshifts is still within the range of obser- 
vational uncertainty. 

Ultimately, we find that our predictions are robust above 
masses Mgai > IO'^Mq, regardless of these changes to our 
model, as the theoretical subhalo mass functions and empir- 
ical halo occupation models are reasonably well-constrained 
in this regime. Below these masses, in any case, it is likely 
that a large fraction of spheroids are relatively small bulges 
in disk-dominated galaxies (of which a large fraction may be 
pseudobulges formed by disk instabilities) and that a large 
fraction of the red galaxy population are satellites (whose red- 
dening may be affected by their mere accretion as a satellite, 
let alone tidal or ram pressure stripping processes, which we 
do not attempt to model). While not dominant in the > 
galaxies with which our modeling is most concerned, these 
processes are certainly important for low-mass populations. 

We further discuss a variety of physical mechanisms that 
may drive the quenching of major merger remnants. In nu- 
merical experiments, the star formation rates of isolated disks 
(i.e. ones cut off from any gas accretion) decay slowly, and 
the galaxies do not move to the red sequence in times < a few 
Gyr (despite allowing for secular instabilities in these sim- 
ulations). This alone is a consideration which should be of 
concern in secular or pure halo quenching models - without 
mergers or some other driver of violence in the system, these 
systems do not efficiently transition to the red sequence in 
the first place. However, it is clear that merger remnants effi- 



ciently exhaust gas and redden rapidly onto the red sequence, 
even without the inclusion of feedback effects (although these 
may be necessary to fully terminate star formation in the most 
high-redshift, gas-rich systems). It is clear that mergers easily 
accomplish the "transition" to the red sequence, even if only 
temporarily. The more difficult question is how such systems 
might prevent future cooling, in order to remain quenched for 
significant periods of time. 

There are, however, a number of feedback sources di- 
rectly associated with major mergers, including purely kine- 
matic "stirring," tidal heating, and shock effects, long-lived 
starburst-driven winds, and (potentially) impulsive, quasar- 
driven outflows. We demonstrate in numerical simulations 
(and from simple scaling arguments) that the combination of 
these feedback effects (even with relatively mild prescriptions 
for their strength) is sufficient to heat several times the initial 
baryon content of the host halo at the time of the merger to 
very high temperatures, at which the cooling time becomes 
longer than a Hubble time. For the ^ 1 /2 of the present ^ 
red galaxy population that has moved onto the red sequence 
since z ~ 1, this single feedback event is sufficient to prevent 
all but < 1% of the galaxy mass from cooling back onto the 
galaxy by z = 0, i.e. sufficient to ensure the galaxy remains 
"red and dead." The problem, however, is potentially more 
severe at high redshifts. Not only must the suppression of 
cooling act for a longer period of time, but a massive halo at 
high redshifts z > 2 may typically grow by a large amount 
(more than an order of magnitude) in mass by z = 0, imply- 
ing that the total baryon content for which cooling must be 
suppressed is larger than that of the galaxy. Moreover, the 
increased densities at these times further suppress the prop- 
agation of feedback-driven shocks, and enhances the cooling 
rates by large factors (^ 100 at z ~ 3 -4). 

We therefore propose a "mixed" solution, in which tradi- 
tional modes of quenching and feedback in quasi-static "hot 
halos" remain the key to suppressing star formation in mas- 
sive systems, but that these are supplemented by mergers, 
which can effectively quench star formation in lower-mass 
systems before these cross the hot halo threshold. The ma- 
jor merger can supplement this traditional quenching mode 
in two ways: first, by temporarily suppressing cooling un- 
til the system naturally develops a hot halo (i.e. crosses the 
quenching mass threshold). This is, even for high redshift 
systems, often much less than a Hubble time (~a couple 
Gyr), and is relatively easily accomplished by the feedback 
effects described above. Second, the strong shocks from the 
merger-driven feedback can accomplish what an accretion 
shock would in a more massive system - i.e. they can create a 
quasi-static hot halo even in low-mass systems which would 
not (independent of a major merger) develop such a halo on 
their own. 

We demonstrate both using simple scaling arguments and 
numerical simulations including feedback, cooling, star for- 
mation, and realistic shock mechanisms, that even conserva- 
tive feedback prescriptions will shock most of the gas within 
the virial radius to temperatures and entropies where the cool- 
ing time becomes much longer than the free fall or dynami- 
cal/compression timescales (the traditional definition of a hot 
halo). Once this hot halo is established inside 7?vir, a quasi- 
static, pressure supported equilibrium is established against 
which newly accreted gas will shock and add to at large radii 
(regardless of the mass subsequently accreted). The energet- 
ics of merger-triggered feedback are sufficient to achieve this 
in all halos Mhaio ^ 1O'^M0 (including halos below the tra- 
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ditional quenching mass threshold), with little dependence on 
redshift (at least from z = 0-6). Once a hot halo is devel- 
oped, the problem of maintaining that hot halo (i.e. prevent- 
ing cooling flows) is no different from the traditional cool- 
ing flow problem (which we are not directly attempting to ad- 
dress here), but the merger remnant already, by definition, has 
the means to heat the halo and supplement it with feedback 
- namely, a relatively massive spheroid and BH which will 
be accreting at low rates (i.e. the ideal seed for "radio-mode" 
feedback). 

This ultimately simple variation on the traditional models 
of quenching in massive systems appears to yield a number 
of qualitatively different predictions, as described above, and 
merits further study. Although, in order to limit the physi- 
cal assumptions being studied, we did not adopt a full semi- 
analytic model, it will be valuable for future studies and com- 
parison to observations to implement such models. There are 
a number of prescriptions one might consider, with varying 
degrees of complexity, which may yield different, testable ob- 
servational predictions. Ideally, such models should consider 
a variety of prescriptions for quenching, and compare the re- 
sults in order to determine what (if any) observational tests 
might break the degeneracies between them. 

(1) Pure Merger Quenching: This is the simplest possi- 
ble model, similar to what we have assumed in this work, 
assuming that a major merger completely suppresses future 
cooling/star formation. Equivalently, one could adopt some 
bulge-to-disk ratio ab ove which cooling is suppressed, as in 
(ICattaneo et al.l l20()6h . or a bulge mass threshold (Mbuige ^ 
3 X IO'^Mq, as in Naab et al. 2006b). 

(2) Merger Feedback/Strong Shocks: Rather than fiat 
quenching, one could allow for some large energy injection 
from feedback (presumably owing to triggered quasar and 
starburst activity) in a merger, and assume that the appro- 
priate shocked quantity of gas has its cooling suppressed, 
or is ejected from the host galaxy and reheated to the halo 
virial temperature (Somerville et al. 2007). This is similar 
to the calculation in Scan napieco & Oh (.2004 *). who demon- 
strate that such an assumption is sufficient to produce down- 
sizing trends below z ~ 2. There are a number of ana- 
lytic models which have been proposed for the effects of 
this feedback, including the blastwa v e mod el calibrated to 
simulations in [H opkins & Hernquisi (l2006l) . the model of 
IScannapieco & O h (2004) in terms of the post-sho ck e ntropy, 
and the temperature/cooling time calculations in § 14.21 

(3) Merger-Induced "Hot Halos": Based on the arguments 



above, it is straightforward to assume that feedback from 
quasar and/or starburst activity triggered in a major merger 
drives the host halo to the quasi-static, hot halo regime. What- 
ever the treatment in the semi-analytic model is for such 
hot halos (i.e. whether they are assumed to be quenched, 
or whether various AGN feedback modes are considered for 
"maintenance" purposes), the host halos of major merger rem- 
nants would be treated identically. 

(4) A "Full Model": Ideally, semi-analytic models could 
incorporate all of the effects above. Based on energetic 
arguments or the simple scaling arguments in § 14.21 or 
adopting some analytic model for a feedback-driven shock 
(IScannapieco & Ohll2004t iHopkins & Hernquisi l2006i) . one 
can calculate the appropriate effects on halo gas. It is then 
possible to consider whether this moves the halo into the hot 
halo regime, or "buys time" until the halo experiences ac- 
cretion shocks and falls into such a regime itself. Feedback 
from low-luminosity AGN, or cyclic accretion inside a hot 
halo, would be allowed, and could further suppress subse- 
quent cooling. 

Future study using high-resolution numerical simulations 
will be essential to ultimately understanding the interplay of 
these complex feedback processes. Simulations with the dy- 
namic range to simultaneously resolve the relevant galactic 
structure and feedback processes and cosmologically rare, 
massive populations are not yet feasible; however, the effects 
of these processes in represent ative systems can be studied in 
detailed zoom-in simulations ( Li et al.ll2006h . Examining, for 
example, the effects of feedback on clumpy accretion at high 
redshift or the details of how merger-driven shocks transform 
the halo cooling structure will be critical to inform theoretical 
models of how these systems quench and suppress cooling 
over cosmic time. The combination of detailed simulations 
used to study the effects of feedback and cosmological models 
which enable predictions for the broad statistical properties of 
rare populations should allow future observations to break the 
degeneracies between different quenching models and tightly 
constrain the history of massive galaxy formation. 
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